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Abstract 



We consider the stochastic heat equation of the following form 
d 

— u t (x) = (Cu t )(x) + b(u t (x)) + a(u t (x))F t (x) for t > 0, x € R d , 
at 

where £ is the generator of a Levy process and F is a spatially-colored, 
temporally white, gaussian noise. We will be concerned mainly with the 
long-term behavior of the mild solution to this stochastic PDE. 

For the most part, we work under the assumptions that the initial data 
uq is a bounded and measurable function and a is nonconstant and Lipschitz 
continuous. In this case, we find conditions under which the preceding 
stochastic PDE admits a unique solution which is also weakly intermittent. 
In addition, we study the same equation in the case that Cu is replaced by 
its massive/dispersive analogue Cu — Xu where A G R. And we describe 
accurately the effect of the parameter A on the intermittence of the solution 
in the case that cr(u) is proportional to u [the "parabolic Anderson model"]. 

Furthermore, we extend our analysis to the case that the initial data 
uq is a measure rather than a function. As it turns out, the stochastic 
PDE in question does not have a mild solution in this case. We circumvent 
this problem by introducing a new concept of a solution that we call a 
temperate solution, and proceed to investigate the existence and uniqueness 
of a temperate solution. We are able to also give partial insight into the 
long-time behavior of the temperate solution when it exists and is unique. 

Finally, we look at the linearized version of our stochastic PDE, that is 
the case when a is identically equal to one [any other constant works also]. 
In this case, we study not only the existence and uniqueness of a solution, 
but also the regularity of the solution when it exists and is unique. 
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CHAPTER 1 



Introduction and Statements of Main Results 

The principle aim of this paper is to describe the asymptotic large-time 
behavior of the mild solution u := {ut(%)}t>o,xefl d °f the stochastic heat 
equation, 

(1.1) ^u t {x) = (Cu t )(x) + b{u t (x)) + a{u t (x))F t (x), 

where t > and x G R d , and the preceding stochastic PDE can be under- 
stood in the sense of Walsh [Wal86] . 

For the most part, we consider the case that the initial data uq is a 
nonrandom, as well as bounded and measurable, function. But we will also 
consider the physically-interesting case that uq is a nonrandom finite Borcl 
measure on R rf . The latter case will be the subject of Chapter 6. 

Throughout we consider only functions a, b : R — > R that are non- 
random and Lipschitz continuous. Also, we let C be the L 2 -generator of a 
d-dimensional Levy process X := {X t } t >o, and assume that X has transition 
functions. 

In the above discussion, we have used the standard notation of probabil- 
ity theory: Namely, gt denotes the evaluation of a [random or nonrandom] 
function g at time t, and never the time derivative of g. This notation will 
be used throughout the rest of the paper. 

As regards the forcing term F in (1.1), we assume that F is a generalized 
Gaussian random field [GV77, Chapter 2, §2.4] whose covariance kernel 
is So(s — t)f(x — y), where the "correlation function" / is a nonnegative 
definite, symmetric, and nonnegative function that is not identically zero. 1 
Alternatively, one can use the following 



The symbol "/" is reserved for this correlation function here and throughout. We never 
refer to any other function as /. 
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in the sense of generalized random fields, where F is a centered generalized 
Gaussian random field with covariance kernel 



(1.3) 



Cov ( / (PdF , / (dF 

'R + xR d JR + xR d 



/ ds f dx [ dy 4> s (x)C, s {y)f(x 
JO JR.- 1 JR d 



v), 



where f <fr dF and f (" dF are Wiener integrals of R + x R d B (s , x) i-> 
4> s (x) and R+ x H d B (s,x) i— >■ £ s (x) with respect to F, and and £ are 
nonnegative measurable functions for which the right-most multiple integral 
in (1.3) is absolutely convergent. 

According to the Bochner-Schwartz theorem [GV77, Theorem 3, p. 
157]: 

(a) The Fourier transform / of / is a [nonnegative Borel] tempered 
measure on R d ; and 

(b) Conversely, every tempered measure / on H d is the Fourier trans- 
form of one such correlation function /. 

The measure / is known as the "spectral measure" of the noise F. Through- 
out, we assume without further mention that F "has a spectral density." 
That is, 

(1.4) / is a measurable function. 

This implies that / is locally integrable on R rf as well. Strictly speaking, 
these conditions are not always needed in our work, but we assume them 
for the sake of simplicity. 

By enlarging the underlying probability space, if need be, we introduce 
an independent copy X* := {X^}t>o of the dual process —X. We can then 
use X* to define a symmetric Levy process X := {X t } t >o on R rf via the 
assignment 

(1.5) X t :=X t + XZ foralH>0. 

Motivated by the works of Kardar, Parisi, and Zhang [KPZ86] and Kardar 
[Kar87], we may refer to X as the replica Levy process corresponding to X 
and will therefore call the resolvent {R a }a>o of X, the replica resolvent. 2 
We will consider the condition that the correlation function / has finite 
a-potential at zero for all a > 0. That is, we consider the following: 



2 These quantities are defined in more detail in Chapter 2. 
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Condition 1.1. (i? Q /)(0) < oo for all a > 

The above condition will imply an existence and uniqueness result for 
the stochastic heat equation (1.1). Moreover, our proof of existence and 
uniqueness is closely linked to the large-time behavior of the solution itself 
[via a priori estimates]. We describe these results next. But first, let us de- 
fine two important quantities: The first denotes the upper LP{P)-Liapounov 
exponent of the solution u := {ut(x)} t>0xe jid to (1.1) at the spatial point 
x e H d : 

(1.6) %{p) := limsup ylnE(|ut(x)| p ) ; 

and the second the upper maximum LP (P)-Liapounov exponent: 

(1.7) 7„(p) := limsup \ sup InE (\u t (x)\ p ) . 

The above two quantities are variants of the well known Liapounov exponent. 
We are now ready to state the first main contribution of this paper. 

Theorem 1.2. Assume that Condition 1.1 holds, and suppose uq : R d — >■ 
R is bounded and measurable. Then, (1.1) has an a.s. -unique mild solution 
which satisfies the following: For all even integers p > 2, 

(1.8) %(p) <inf{/3>0: Q(p,0)<l}, 
where 

pLip b 



(1-9) Q(p,P) ■= ^f± + z p Up a ^(R 2f3/p f)(0), 

and z p denotes the largest positive zero of the Hermite polynomial He p . 

Let us make two remarks before we continue with our presentation of 
the main results of this paper. The first one is consequence of the above 
result. 

Remark 1.3. It is possible to deduce from Condition 1.1 and the mono- 
tone convergence theorem that lim a _ s . oo (^ a /)(0) = 0. This, in turn, implies 
that 

(1.10) lim Q(p,0) = O. 

/3-s-oo 

Consequently, Theorem 1.2 implies among other things that 7„(p) < oo for 
ah>G(0,oo). □ 
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Remark 1.4 (Borrowed from [FK09, Remark 2.2]). It might help to 
recall that 

(1.11) Hek(x) = -^j^Hk (^^j f° r an i nte § ers ^ > and x G R, 
where {Hk}^L is defined uniquely via the following: 

OO 

(1.12) e -2^-t 2 = V TT #fc(x) for alU > and x G R. 

rCI 

fc=0 

It is not hard to verify that 



(1.13) «a = 1 and z 4 = \J 3 + y/E « 2.334. 

This is valid simply because He2(x) = x 2 — 1 and iie^x) = x 4 — 6x 2 + 3. In 
addition, 

(1.14) z p ~ as p — )■ oo, and sup ( — — ) = 2; 

P >1 VVP/ 

see Carlen and Kree [CK91, Appendix]. □ 



Next we put Theorem 1.2 in the context of the existing literature on the 
stochastic heat equation. With this in mind, define for all j3 > 0, 



(1.15) T ^ : =^/" 



/(0 



/3 + 2Re*(0 

where * is the characteristic exponent of the Levy process X. Dalang 
[Dal99] has established a very general result which guarantees the existence 
and uniqueness of solutions to large families of SPDEs. If we apply Dalang's 
result to the present parabolic problem (1.1), then we find the following: If 
uq is a constant, then the condition 

(1.16) T(l) < oo 

insures the existence and uniqueness of a [mild] solution to (1.1). Moreover, 
Dalang's result shows that (1.16) is necessary and sufficient for existence and 
uniqueness in the case that (1.1) is a linear SPDE; that is, when cr(u) = 1 and 
b(u) = 0. For closely-related results [that also include hyperbolic equations] 
see Carmona and Molchanov [CM94], Conus and Dalang [CD08], Dalang 
and Frangos [DF98], Dalang and Mueller [DM03], Dalang, Mueller, and 
Tribe [DMT08], Dalang and Sanz-Sole [DSS09], and Peszat and Zabczyk 
[PZ00]. 
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Our next result implies among many other things that Dalang's con- 
dition (1.16) is generically equivalent to the potential-theoretic Condition 
1.1. 

Theorem 1.5 (A Maximum Principle). For all (3 > 0, 

(1.17) (RfsfM = sup (Rfsf)(x) = T(/3). 

x€R d 

Thus, Condition 1.1 holds if and only if (1.16) holds. Furthermore, if Con- 
dition 1.1 [and/or (1.16)] holds and f is lower semicontinuous, then for all 
(3 > there exists irp G Co(R d ) such that Rpf = Tip almost everywhere. 

In light of Theorem 1.5 and Dalang's theorem [Dal99], the novel con- 
tributions of our Theorem 1.2 are: 

(a) The condition that uo is a constant can be improved to one about 
the boundedness of uq [this can also be derived by adapting the 
method of Dalang [Dal99] to the present setting]; and more signif- 
icantly 

(b) We obtain a uniform upper bound for the maximum L p (P)-moment 
Liapounov exponent of the solution to (1.1) as an a priori conse- 
quence of the existence of the solution. 

This second contribution leads to the weak intermittence of solutions, which 
is a notion that is rooted in the literature of statistical mechanics. With this 
in mind, let us recall the following [FK09]: 

Definition 1.6. Suppose that there exists an a. s. -unique solution u := 
{«((i)}( >0a , eR d to (1.1). We say that u is weakly intermittent if 

(1.18) 7*{p) < oo for all p G [2,oo), and inf J x (2) > 0. 

xeR d 

The same reasoning that was employed in [FK09] can be used to deduce 
that if the solution to (1.1) is nonnegative for all t > 0, then weak inter- 
mittence implies the much better-known property of intermittency [CM94, 
Mol91, ZRS90]; that is, the property that 

T \P) rl 

(1.19) p H> — is strictly increasing on [2 , oo) for all x G R . 

There is a large literature which shows that, under further mild hypotheses 
on C and/or /, if no is nonnegative then the solution to (1.1) is nonnegative 
at all times; see, for example the papers by Assing and Manthey [AM95], 
Carmona and Molchanov [CM94], Donati-Martin and Pardoux [DMP93], 
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Hausmann and Pardoux [HP89], Kotelenez [Kot92], Manthey [Man86], 
Manthey and Stiewe [MS92, MS91], Mueller [Mue91], Nualart and Par- 
doux [NP92], and Shiga [Shi97, Shi94]. 3 

Thus, we can draw the conclusion that, in all such cases, weak intermit- 
tence actually implies intermittency. 

A quick calculation, using only Holder's inequality, shows that p i— >■ 
J x (p)/p is always nondecreasing on [2,oo). However, the mentioned strict 
monotoncity does not always hold. When it does hold, then it has some 
physical significance; see Zeldovitch, Ruzmaikin, and Sokoloff [ZRS90] for 
a physical discussion of intermittency. And Molchanov [Mol91] for a math- 
ematical explanation of that physical phenomenon. 

Our next main goal is to find nontrivial conditions that guarantee the 
weak intermittence of the solution to (1.1). In light of Theorem 1.2, we 
aim to derive a positive lower bound on inf^gRd 7 X (2). Unfortunately, it is 
quite hard to do this at the level of generality of the conditions of Theorem 
1.2. In fact, informal arguments suggest that the solution to (1.1) might not 
always be weakly intermittent. Thus, we seek to find reasonable restrictions 
of the various parameters of (1.1) which guarantee that the solution to (1.1) 
is weakly intermittent. 

Let g denotes the Fourier transform of a locally-integrable function g, 
and consider the following: 

Condition 1.7. Suppose: 

(1) f{0 depends on £ G R d only through |£i|, . . . , \£ d \; 

(2) \£j\ i->- /(£) is nonincreasing for every j = 1, . . . , d; and 

(3) Re#(£) depends on £ G R d only through |£i|, . . . , 

These are relatively mild provisions on the spectral density / and the 
process X. Our conditions on the spectral density can be applied to all of 
the examples that we would like to cover. It is possible to show that they 
include the following choices for /: 

(i) Ornstein-Uhlenbeck-type kernels. 
f(x) = Cl e- C2Ma 



3 In connection to matters of positivity and regularity, we mention also a closely-related 
and fundamental paper by Dawson, Iscoe, and Perkins [DIP89], where (1.1) with cr(w) = 
const • ^Ju is considered. And positivity of the solution is shown to follow from many- 
particle approximations to the underlying SPDE. 
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for constants c±, c 2 G (0 , oo) and a G (0,2]; 
(ii) Poisson kernels. 



fix) 



(IMI 2 +C 2 ) (<m)/2 

for ci, C2 G (0 , oo). 
(iii) Cauchy kernels. 

ci 



/(£ ) = const • e 



-const 



n-=i(c 2 +x]) 



= const . e - const ^=i fcl 



for ci, C2 G (0 , oo); and 
(iv) ffiesz kernels. 

f^) = 7^77 



/(0 



const 



for c G (0 , oo) and a G (0 , <f). 

And one can construct a great number of other permissible examples as well. 

Having introduced Condition 1.7, we can now present the third main 
result of this paper. 

Theorem 1.8. Suppose 6 = and Conditions 1.1 and 1.7 hold. Suppose, 
in addition, that rj := inf xeR d Uoix) > and there exists L a G (0 , oo) such 
that aiz) > L a \z\ for all z G R. Then, 



(1.20) 



inf 7 X (2) > sup /3 > : (^/)(0) > 
where sup := 0. 



Before we pause to make a few remarks, let us briefly study an example. 
Consider the case that (-Ro/)(0) = oo. In that case, (Rpf)(Q) > 2 d ~ 1 /h 2 a for 
all /3 > sufficiently small. Hence, in this case, the hypotheses of Theorem 
1.8 guarantee weak intermittence of the solution to (1.1) without further 
restrictions. 

Remark 1.9. (1) In the case that F is space-time white noise, the 
condition "cx(,z) > L^z)" can be replaced with the slightly-better 
condition a \<r(z)\ > K\z\ n [FK09]. 
(2) We will see later on that, when d = 1, the lower bound (1.20) 
and the upper bound (1.8) can sometimes match. However, the 
two bounds can never agree when d > 2. This phenomenon is due 
to the fact that level sets of /? >->■ (^/)(0) cannot describe the 
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growth of u exactly. The correct gauge appears to be a much more 
complicated function, except in the cases that F is space-time white 
noise and when d = 1; compare with [FK09] for results on the case 
that F denotes space-time white noise. □ 

We are aware of a few variants of Theorem 1.8, but the next one is 
perhaps the most striking since it assumes only that the nonlinearity term 
a is asymptotically sublinear. Thus, the local behavior of a is shown to not 
have an effect on weak intermittence, provided that the initial data uq is 
sufficiently large. A significant drawback of this result is that its proof does 
not provide any information about how large "sufficiently large" should be. 
We introduce the following condition. 

Condition 1.10. (-R o /)(0) = oo. 

The following result is the mentioned variant of the Theorem 1.8. 

Theorem 1.11. Suppose 6 = and Conditions 1.1, 1.7, and 1.10 hold. 
Suppose, in addition, that a > pointwise, and q := liminfui^^ cr(z)/|z| > 
0. Ifu (x) > and P{u t (x) > 0} = 1 for all t > and x £ K d , then 

(1.21) 73,(2) > for every x G K d , 
provided that n := inf^gj^d uq(x) is sufficiently large. 

The preceding results describe our main contributions to the analysis 
of the stochastic heat equation (1.1) in the case that a is not a constant 
and that uq is a bounded and measurable function. But we also study the 
linearization of (1.1); this is the case when a is identically equal to one. In 
addition to studying existence-and-uniqueness issues, we use the theory of 
Gaussian processes to study continuity properties of the solutions. Moreover, 
we produce a class of interesting examples which we briefly describe next. 

Consider the linear stochastic heat equation 

(1.22) - Ut (x) = (Au t )(x) + F t (x), 

where uq = 0, x G R 3 , t > 0, and the Laplacian acts on the x variable 
only. Then, we construct families of noises F which ensure that (1.22) has a 
solution u := {ut(^)}t>o,a;GR d that is a square-integrable random field. But 
that random field is discontinuous densely with probability one. In fact, 
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outside of a single null set [of realizations of the process u] , 

(1.23) sup ut{x) = — inf u t (x) = oo, 

(t,x)ev (t,x)ev 

for all open balls V C R+ x ~R d with rational centers and radii! 4 We 
know of only a few examples of SPDEs with well-defined random-field solu- 
tions that have unbounded oscillations densely; see Dalang and Leveque 
[DL06, DL04b, DL04a], Mytnik and Perkins [MP03], and Foondun, 
Khoshnevisan, and Nualart [FKN09]. The preceding (1.23) yields a quite- 
simple example of an otherwise physically-natural stochastic PDE [the op- 
erator is the Laplacian in R 3 and the noise is white in time] which has a 
very badly-behaved solution. 

This paper was influenced greatly by the theoretical physics literature 
on the "parabolic Anderson model" (see, for example, Kardar, Parisi, and 
Zhang [KPZ86], Krug and Spohn [KS91, §5], Medina, Hwa, Kardar, and 
Zhang [MHKZ89], and the book by Zeldovitch, Ruzmaikin, and Sokoloff 
[ZRS90]), as well as the mathematical physics literature on the very same 
topic (see, for example, Bertini and Cancrini [BC98, BC95], Bertini, Can- 
crini, and Jona-Lasinio [BCJL94], Bertini and Giacomin [BG99, BG97], 
Carmona, Koralov, and Molchanov [CKM01], Carmona and Molchanov 
[CM95, CM94], Carmona and Viens [CV98], Cranston and Molchanov 
[CM07b, CM07a], Cranston, Mountford, and Shiga [CMS05, CMS02], 
Florescu and Viens [FV06], Gartner and den Hollander [GdH06], Gartner 
and Konig [GK05], Hofsted, Konig, and Morters [vdHKM06], Konig, La- 
coin, and Morters [KLMS09], Lieb and Liniger [LL63], Molchanov [Mol91], 
and Woyczyhski [Woy98] for a partial listing). Furthermore, there are in- 
teresting variations of the parabolic Anderson model that correspond to 
continuous directed-polymer measures; see Comets and Yoshida [CY05] 
and Comets, Shiga, and Yoshida [CSY04, CSY03]. 

In a nutshell, the parabolic Anderson model is equation (1.1) where 
a(u) is proportional to u. There are many good reasons why that equation 
has been studied intensively; see for instance the Introduction of Carmona 
and Molchanov [CM94]. Two such reasons are that the parabolic Anderson 
model it is exactly solvable in the two cases where uq = constant and uq = Sq; 



It would be wonderful to construct versions of these examples that apply to fully nonlinear 
problems; but we do not know how to proceed in a fully nonlinear [or even semilinear] 
setting. 
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and it is related deeply to the stochastic Burgers equation as well as the KPZ 
equation of statistical mechanics. 

And perhaps not surprisingly, the results of our Theorems 1.2, 1.8, and 
1.11 are sharpest for the parabolic Anderson model, particularly when d = 1. 
However, an inspection of Theorems 1.2 and 1.8 reveals an inconsistency: 
Our upper bound on the Liapounov exponent [Theorem 1.2] does not require 
the drift b to be zero; whereas our lower bound [Theorem 1.8] does. 

David Nualart has asked us whether we know how the drift b can affect 
the weak intermittence of the solution to (1.1). This seems to be a hard ques- 
tion to answer rigorously when the drift 6 is a general Lipschitz-continuous 
function. But it is intuitively clear that a sufficiently-strong drift ought to 
destroy the natural tendency of the solution to be weakly intermittent. 

Although we are not aware of general theorems of this type, we are able 
to give a partial answer to Nualart's question; and the striking nature of 
that partial answer confirms our initial suspicion that it might be rather 
difficult to answer D. Nualart's question in good generality. 

Here is an instance where we can rigorously prove weak intermittency: 
Consider the one-dimensional parabolic Anderson model for the relativistic 
[or massive/dissipative] Laplacian; i.e., the stochastic PDE 

d X 

(1.24) inMx) = {Au t )(x) + - u t (x) + Ku t (x)F t {x), 
Ot 2 

where t > and x G R, k ^ 0, A G R, and uq : R — > R is a measurable 

function that is bounded uniformly away from zero and infinity. Let us 

consider the special case that the correlation function of the noise is of 

Riesz type; that is, 

(1.25) f(z) := ||z||- 1+b for all z G R, 

where b G (0,1). Then, Example 5.8 on page 71 implies that weak intermit- 
tence holds if and only if 

T(6/2)r((b + l)/2) l 2/(1+6) 

This completes our investigation of the solution to (1.1) when the initial 
data no is a bounded and measurable function. 

We conclude this paper by considering (1.1) in the other physically- 
interesting cases where no is a finite Borel measure on R rf . This condition 
on no is very natural, as in many physical applications no denotes the ini- 
tial distribution of particles in a particle system in a disordered medium 



(1.26) A>-| K | 4 /( 1+6 )8-( 1 - 6 )/( 1+b ) 
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[KPZ86, MHKZ89, KS91, Mol91, CM94, BC95]. From a purely 
mathematical point of view, this problem is interesting since it leads us 
to a different notion of a solution, which we call temperate. Our proposed 
temperate solutions differ from the much better-known notion of "mild solu- 
tions" [Wal86, Chapter 3]. They also lead us to an extension of the notion 
of a stochastic convolution, which defines stochastic integrals at almost ev- 
ery time, rather than pointwise. We believe that the notion of temperate 
solutions has other uses in describing otherwise hard-to-define SPDEs. In 
order to describe our results on temperate solutions — the final main con- 
tributions of this paper — we would have to develop some machinery. The 
details can be found in Chapter 6. 

A brief outline of the paper follows: In Chapter 2 we review, very briefly, 
some analytical facts about Levy processes and their generators, and also 
construct examples that will be used in subsequent chapters. 

Chapter 3 is concerned with positive-definite functions and their connec- 
tions to potential theory and harmonic analysis. And Theorem 1.5 is shown 
to be a consequence of these connections. Chapter 3 also contains a prob- 
abilistic characterization of the analytic condition (1.10) and Condition 1.1 
in terms of continuous additive functionals of the replica process X. Also, 
a family of useful correlation functions is constructed in that chapter; that 
construction uses the results of Chapter 2 on probabilistic potential theory. 

In Chapter 4 we study the linearization of (1.1), and derive necessary and 
sufficient conditions for the existence and spatial continuity of the solution. 
We also consider various examples that include (1.22) above. 

Chapter 5 contains the proofs of Theorems 1.2 and 1.8. In that chap- 
ter we consider also the relativistic version of (1.1), thereby constructing 
examples that include the mentioned analysis of (1.24). 

Finally, we conclude with Chapter 6, where a new notion of "temperate 
solution" is introduced. And then that notion is used to produce solutions 
to (1.1) in the case that 6 = and uq belongs to a suitable family of finite 
Borel measures on RA 

Let us conclude the present chapter by introducing some notation that 
will be used throughout the paper. For all integers k > 1, 

(1.27) ||x|| := {x\ -\ h x 2 k ) 1/2 for every x G R fc . 



12 1. INTRODUCTION AND STATEMENTS OF MAIN RESULTS 

And if g : R k -»• R is a function, then 

(1.28) Lip 9 := sup \9(x)-9(y)\ _ 

x,yeR. k W X y\\ 

This socalled Lipschitz constant of g is well defined, but might be infinity. 

Throughout this paper, " ~ " denotes the Fourier transform in the sense 
of L. Schwartz; our Fourier transform is normalized so that 

(1.29) 9(0 ■■= I e ix <g(x)dx for all £ G R rf and g G L 1 (R (i ). 

Finally, if g : K d R is a function, then we define 

(1.30) g(x) := g(-x) for all x G R d . 

And similarly, if jj, is a Borel measure on R d , then for all Borel sets A C R d , 

(1.31) /t(A) := n(-A), where -^4 := {-a : a£ A}. 
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CHAPTER 2 



Levy Processes 

2.1. Preliminaries 

We begin this chapter with the definition of a Levy process which through- 
out this paper, will be denoted by X := {X t }t>o- 

Definition 2.1. We say that X := {Xt}t>o is a Levy process if: 

(1) X t + S — X s is independent of the sigma-algebra generated by {^r}>e[o,.s] 
for every s, t > 0; 

(2) X t + S — X s has the same distribution as X t for every s, t > 0; 

(3) t ^ X t is continuous in probability; that is, X s converges to X t in 
probability as s — > t; and 

(4) X = 0. 

By adopting a suitable modification of the paths, we can and will always 
assume, without loss of generality, that the trajectories of X are cadlag; i.e., 
t i->- X t is almost surely right-continuous with left limits. Comprehensive 
treatments can be found in the books by Bertoin [Ber96], Jacob [Jac05], 
Kyprianou [Kyp06], and Sato [Sat99]. 

Let vrit denote the distribution of X t for every t > 0; that is, 

(2.1) m t {A) := P{X t G A} for all t > and Borel sets A C R d . 

Let us recall that throughout this paper we are assuming that the process 
X has transition functions; that is, 

(2.2) m t (dx) < dx for all t > 0. 

According to Theorem 2.2 of Hawkes [Haw79], we can always select a ver- 
sion of these transition functions that has the following regularity features: 

(1) j A Pt(z) dz = m t (A) for all t > and Borel sets A C R d ; 

(2) (0 , oo) x R d 9 (t , x) i-)- pt(x) G R + is Borel measurable; 

(3) x i->- pt(x) is lower semicontinuous for all t > 0; 

(4) p t+s (x) = (pt *p s ){x) for all s,t > and x G R d , 
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where "*" denotes the convolution operator, defined in the sense of L. 
Schwartz. We work only with such a version of these transition functions. 
Note that for all t > and x G Ti d , and for every Borel- measurable function 
4> : K d -»• R+, 

E(j)(x + X t )= / <f>(z)pt (z - x) dz 

(2.3) Jr* 

= {4>*pt){x), 

where we recall pt(x) := pt(—x). 

Alternatively, one can work with the semigroup {Pt}t>o of X, which is 
defined via 

(2.4) {P t <t>){x):=V<t>{x + X t ). 

It is easy to verify that {Pt}t>o is a Feller semigroup; i.e., 

(2.5) P t : C (R d ) -> C (R d ), 

where Co(R rf ) denotes the collection of all continuous functions g : H d — > R 
that vanish at infinity. In fact, under the present conditions, {Pt}t>o is 
strong Feller in the sense of Girsanov [Gir60]; see Hawkes [Haw79]. 
Let us emphasize that 

(2.6) P t (t> = <j)*p u 

valid for all t > 0, and for instance for every nonnegative Borel-measurable 
functions (j) : H d — > R. 

Let {R a } a >o denote the resolvent of {.Pt}t>o; i-e., 

POD 

(2.7) R a := / e~ as P s ds. 

Jo 

It follows that if <f) : R rf — > R+ is Borel measurable, then 
(2.8) 



poo 

(R a (/))(x) = / 4>(z)r a (x - z) dz 
Jo 



= (4>*r a )(x), 

where 

POO 

(2.9) r a (x) := / e- at p t (x) dt for a > and x G TL d . 

Jo 

Each "a-potential density" r a (x) is well defined, but could well be infinity 
at some [in fact, even all, when a = 0] x € H d . Nevertheless, the reg- 
ularity properties of the transition functions imply that every r a is lower 
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semicontinuous. Furthermore, 

(2.10) R a : C (R d ) -> C (R d ) for every q > 0. 

In fact, R a (Co(R d )) is uniformly dense in Co(R d ) when a > 0; see Blumen- 
thal and Getoor [BG68, Exercise (9.13), p. 51]. 

The characteristic exponent of the process X is a function \f : R ' — > C 
that is defined uniquely via 

(2.11) Ee^ Xt = e - **® for all £ € R rf and i > 0. 

The Levy-Khintchine formula [Ber96, Theorem 1.2, p. 13], and a theorem 
of Schoenberg [Sch38a, Sch38b] together imply that the family of all Levy 
processes is in one-to-one correspondence with the family of all "negative- 
definite functions." 



2.2. The Generator 

We will be working with the L 2 -theory of generators, as developed, for 
instance, in the book by Fukushima, Oshima, and Takeda [FOT94] for 
more general Markov processes. We outline the details in the present special 
case; matters are greatly simplified and in some cases generalized because 
of harmonic analysis. 

Define 

(2.12) Dom[£] = {</> G L 2 (R d ) : G L 2 (R d )| . 

Plancherel's theorem guarantees that <fi G Dom[£] if and only if (j) : R rf — > R 
is Borel-measurable, locally integrable, and 

(2.13) / (1 + |^(0| 2 ) |0(O| 2 d^<oo. 
It is well known that the following holds: 

(2.14) limsup^il < oo. 

neiKoo ll£ll 

This can be derived directly from the Levy-Khintchine formula; see the book 
by Bochner [Boc55, (3.4.14), p. 67]. In the "Notes and References" section 
for Chapter 3 [Boc55, p. 169] of his influential book, Bochner ascribes (2.14) 
in part to Kolmogorov and Levy. 
Recall that 

(2.15) W 1,2 (R d ) := G L 2 (R d ) : V</> G L 2 (R d )} . 
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Because of Plancherel's theorem, 

|2 , IIY7JJI2 



(2.16) 

Therefore, we can see from (2.14) that 



l 2 (R d ) + \\v<t>\\h m = / Rd (i + IKII 2 ) \H0\ 2 df. 



(2.17) 



W L > 2 (R d ) C Dom[£] C L (R ). 



Of course, 5 is dense in W 1,2 (R d ), when the latter is endowed with the usual 
Sobolev norm, 



(2.18) \\Hwi,2 {Rd) := |H|£ 2(Rd) + ||V0||£ 2(Rd) . 

According to Plancherel's theorem, 



(2.19) 



1 



" (2vr) d 7 Rd 
for all i > and ip, <p G L 2 (R d ). Moreover, 

' J L 2 (R d ) 



(2.20) 



-t*(0 



d£. 



It follows easily from this and (2.12) that 



(2.21) 



£<f> := lim • 
40 



exists in L 2 (H d ) if and only if € Dom[£]. Indeed, the sufficiency follows 
from the elementary bound 



(2.22) 



e -**fc) _ 1 



<tmo\, 



and the Cauchy-Schwarz inequality. And the necessity follows from Fatou's 
lemma, upon setting tp := (Pt(/> — 4>)/t- 

Thus, we have the socalled generator [L 2 -generator, in fact] C, defined 
on its domain Dom[£]. In addition, C can be thought of as a convolution 
[or pseudo-differential] operator with multiplier [or symbol] C = — VP. More 
precisely, 

(2.23) £0(0 = -* (00(0 for a11 <t> G Dom[£] and £ G R rf . 

Let us note that for all i > 0, £ G R d , and G ^(R 4 ), 

2 



(2.24) 



p*<K0 



-2tRe*(£) 



10(01 ■ 
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Therefore, the well-known nonnegativity of Re^^) — which we prove at the 
beginning of the following section — implies the following. 

Lemma 2.2. Pt is a contraction on W 1>2 (R d ) for all t > 0. Hence, aR a 
is also a contraction on W 1,2 (R, d ) for all a > 0. 



2.3. The Replica Semigroup and Associated Sobolev Spaces 

Let X* denote an independent copy of the Levy process —X and, fol- 
lowing Levy [Lev37], define 

(2.25) X t := X t + XI for all t > 0. 

It is easy to see that X* := {X^} t >o is the dual process to X, and X := 
{X t }t>o is a symmetric Levy process on H d . And if we denote the distribu- 
tion of X t by fht, then 

(2.26) fht(A) = (m t *rh t )(A) for all Borel sets ACTL d , 
where fht (A) := mt(—A). Note that the Fourier transform of fht is 

= \M0\ 2 

( 2 - 27 ) =e -2tRe*«)_ 

Among other things, this implies the classical fact that 

(2.28) Re* (0 > for all £ G K d . 

The absolute-continuity condition (2.2) implies that every fht is abso- 
lutely continuous with respect to the Lebesgue measure on H d [t > 0]. We 
denote the resulting transition density by pt ■ Every p t is a symmetric func- 
tion on K d [t > 0]. 

We can always choose a version of p that has good regularity features 
[of the type mentioned earlier for p\. In fact, the following version works: 

p t (x) := {pt*Pt){x) 

(2.29) r 

= / p t (x + z)pt(z) dz for x G R and i > 0. 

Equivalently, if P := {Pt}t>o denotes the semigroup of X, then 

(2.30) P t = P t Pt for all t > 0, 

where P t * denotes the adjoint of Pt in L 2 (R d ). Every P t is a self-adjoint 
contraction on L 2 (R d ). 
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Motivated by the work of Kardar [Kar87], we refer to X and P respec- 
tively as the replica process and the replica semigroup. The corresponding 
generator is denoted by C and its domain by Dom[£]. 

For all a > 0, wg Ccin define the replica ot-potenticil defisity v a a,s 

POO 

(2.31) f a (x) := e- as p s (x)ds for all x G Ti d . 

Jo 

Clearly, f a {x) is well defined; but f a (x) can be infinite for some [and even 
all, in the case that a = 0] x G RA The resolvent R := {R a }a>o of the 
semigroup P can also be defined as follows 

poo 

(R a <f>)(x) := / e- as (P s <t>)(x)ds 

(2.32) Jo 

= / H z ) f a( z ~ x )dz, 

for all a > and x G R d . Since f a is a symmetric function on H d , it follows 
that R a 4> = <j) * f a . 

The preceding quantity which is called the a-potential of <p makes sense, 
for example, if <j) '■ R d ~^ R+ is Borel measurable, or when (j) G L p (R d ) for 
some p G [1 , oo] because every P s is a contraction on L p (R d ). 



2.4. On the Heat Equation and Transition Functions 



We begin by recalling some generally-known facts about the fundamental 
[weak] solution to the heat [or Kolmogorov] equation for £: We seek to find 
a function H such that for all i > and x G R rf , 

(2.33) m H ' (x) = (CHt)ix) - 

Ho = 5 Z , 

where z G R d is fixed. By rewriting the above in terms of the [spatial] 
Fourier transform H t and using the fact that the [weak] Fourier transform 
of CH t is —ty-H t [see (2.23)], we obtain the ordinary differential equation, 



(2.34) 



Ho(0 



— J£-z 



The unique solution to this ODE is 
(2.35) H t (0 
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Direct inspection of the Fourier transform reveals that Ht(x) = pt(z — x). 
Thus, we find that the fundamental solution to (2.33) is the measurable 
function (0 , oo) x R d x H d 3 (t ; x , y) i-> p t (y — x). In particular, we might 
observe that in order to have a function solution to (2.33), it is necessary as 
well as sufficient that the underlying Levy process X has transition densities. 

We are thus led to the natural question: "What are the necessary and 
sufficient conditions on the characteristic exponent ^ that ensure the exis- 
tence of transition densities of the corresponding Levy processes"? Unfortu- 
nately, there is no satisfactory known answer to this question at this time, 
though several attempts have been made in the first half of the twentieth 
century; see, for example Blum and Rosenblatt [BR59], Fisz and Varadar- 
jan [FV63], Hartman and Wintner [HW42], and Tucker [Tuc65, Tuc64, 
Tuc62]. 

More recently, Bass and Cranston [BC86] applied Malliavin calculus to a 
family of stochastic differential equations driven by jump noises. Their result 
can be used to supply good sufficient conditions that ensure the existence 
of smooth transition functions. And in [NS06], Nourdin and Simon have 
proved, among other things, that if a Levy process has transition functions, 
then so does the same process plus a drift. 

We will use the following unpublished result of Hawkes. It typically 
provides a good-enough sufficient condition for the existence of transition 
functions. We include a proof in order to document this interesting fact. 

PROPOSITION 2.3 (Hawkes [Haw84]). The following conditions are equiv- 
alent: 

(1) Condition (2.2) holds and p t G L 2 (R d ) for all t > 0; 

(2) Condition (2.2) holds and p t G L°°(R d ) for all t > 0; 

(3) Condition (2.2) holds and pt G L 2 (R d ) for almost every t > 0; 

(4) Condition (2.2) holds and pt G L°°(R d ) for almost every t > 0; 

(5) exp(-Re^) G L<(R d ) for all t > 0. 

(6) exp(— Re^) G L'(R d ) for almost every t > 0. 

Moreover, any one of these conditions implies that: (i) (t,x) >->■ p t (x) has a 
continuous version which is uniformly continuous for all (t , x) G [q , oo) x R d 
for every rj > 0; and (ii) p t vanishes at infinity for all t > 0. 
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Proof (Hawkes [Haw84]). Recall that J Rd p t (x) dx = 1 and p t = 
Pt/2 *Pt/2- Therefore, two applications of Young's inequality yield 

(2.36) INIl-(r<*) < lbt/2|li2 (Rd) < ||Pt/2||i,oo(Rd) for all t > 0. 

Consequently, (1)^(2) and (3)<^(4). 

Next let us suppose that (6) holds. Because 

, , m)\ = 

(2.37) 

Plancherel's theorem ensures that 

g \\Pt\\h ( n d) = j^d\\Pt\\h {Rd) 

<- || -2tRe*|| 

- Il e llL 1 (R d ) • 

Since Re^ > 0, it follows from (6) that p t £ L 2 (R d ) for every t > 0; i.e., 
(6)=>(1). Moreover, we have — in this case — the following inversion formula: 
For almost all x € R d and every t > 0, 

(2.39, Ptix) = -±- d J^-«-,-.no i( . 

It remains to prove that (1) and equivalently (2) together imply (5). 
Recall that pt(x) := pt{—x) and observe that 



(2.40) pt /4 * p t/A 



c 



-(t/2)Re* 



Therefore, by Plancherel's theorem, 



I -Re* II 

r llL*(R d ) 



e -(t/2)Re* 

L 2 (R d ) 
\d\\„ .,. si II 2 



( 2 - 41 ) =(2vr) d ||p t/4 *^ /4 ||^ (Rd) 

< ( 2 7r) d |bt/4*Pt/4|| L oo (Rd) - 
Consequently, Young's inequality, implies that 

(2.42) l|e- Re *IL W < (2vr) d ||^/ 4 |li 2(Rd) , 

which has the desired effect. 

Finally, if (5) holds then the inversion theorem applies and tells us that 
we can always choose a version of p that satisfies the properties of the final 
paragraph in the statement of theorem. □ 
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There is also the following 1970 theorem of J. Zabczyk, which charac- 
terizes (2.2) in special though important cases. 

PROPOSITION 2.4 (Zabczyk [Zab70, Example (4.6)]). If d>2 and * is 

a radial function, then (2.2) holds if and only if 

(2.43) lim tt(0 = oo. 

||£||-HX> 

We now put this beautiful result in context. By the Riemann-Lebesgue 
lemma, if the process X has transition functions then equality (2.43) holds. 
Proposition 2.4 states that the converse also holds, provided that d > 2 and 
^ is radial. In other words, the Riemann-Lebesgue lemma is a necessary 
and sufficient condition for the existence of transition functions whenever 
d > 2 and ^ is radial. As was mentioned in Zabczyk [Zab70] , the preceding 
is not in general true when d = 1. This can be seen by considering to 
correspond to two independent one-dimensional Poisson processes. 

2.5. On a Family of Isotropic Levy Processes 

The main result of this section will be needed to construct a counterex- 
ample in Chapter 4. It is possible that it is known but we were not able 
to find an explicit reference. So we provide a complete proof. We begin by 
recalling a few definitions used to study Levy processes. 

We say that a Levy process X := {Xt}t>o is isotropic if its characteristic 
exponent ^ is a radial function [and hence also real- valued and nonnegative] . 
Such processes are also known as radial processes; see Millar [Mil73]. 

A [standard] subordinator r := {rt}t>o is a one-dimensional Levy pro- 
cess that is nondecreasing and To := 0. According to the Levy-Khintchine 
formula [Ber99, Theorem 1.2, p. 13], every subordinator r is determined 
by the formula 

(2.44) Ee" ATt =e"'* (A) , 
where t, A > 0, and 

(2.45) $(A) = J (l- e" Az ) n(de), 
for a Borel measure LT on (0 , oo) that satisfies 

(2.46) / (1 A x) n(dx) < oo. 

The function <I> is the socalled Laplace exponent of the subordinator r. 
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We have the following lemma. 

Lemma 2.5. Choose and fix two numbers p £ (0,1) and q € R. Then, 
there exists a subordinator r on R + whose Laplace exponent satisfies 

(2.47) < inf - - v / < sup - - v < oo. 

v ; A>c AP(logA)?/ 2 A >c AP(logA)9/ 2 



Proof. Define a measure II via 

x -i-p (log(l/x)) 9/2 if < x < ±, 
otherwise. 



, , \ n(dx 

(2.48) 



dx 

Since p £ (0,1), it follows that (2.46) holds, whence II is a Levy measure. 

We can also apply the definition (2.45) of the Laplace exponent and 
write $(A) = X P Q(X), where 

/■A/2 i — x 

(2.49) Q(A) = / —— (\og(\/x)) q / 2 dx for A > 0. 

Jo x 

In order to complete the proof, we will verify that Q(A) x (log \) q / 2 for 
A > e. 1 We do so in the special case that q > 0; similar arguments can be 
used to estimate Q(X) in the case that q < 0. 
Whenever A > e, we can write 

(2.50) Q(X) := h + I 2 , 
where 

PA/2 j 



(2.51) 

Evidently, 
(2.52) 



/"A/2 i _ -X 

h: = / - T T^r(log(A/x)) ( '/ 2 dx, 



/c 



dx 



X i+ P 

-(logA)«/ 2 . 

Since I\ > 0, it remains to prove that 

(2.53) / 2 ?<(logA) 9/2 forA>l. 



^As usual, h(x) x g(x) over a certain range of x's is short-hand for the statement that, 
uniformly over that range of x's, h(x)/g(x) is bounded above and below by positive and 
finite constants. 



2.5. ON A FAMILY OF ISOTROPIC LEVY PROCESSES 



23 



We establish this by deriving first an upper, and then a lower, bound for 1-2,. 
Because 1 — exp(— y) < y for y > 0, and since sup 2e ( ,i) z e log(l/z) < oo for 
all f £ (0,1), it follows that 



1 / 1™71 /^\«/ 2 



( ^ log(l/x) V 7 



(2.54) " ./„ V. log A J X p 

< const • (log\) q/2 . 

And a similar lower bound is obtained via the bounds: (i) 1— exp(— x) > x/2; 
and (ii) log(A/x) > log A; both valid for all x£ (0,1). □ 

The following is the main result of this section. It gives a special con- 
struction of an isotropic Levy process X := {X t } t >Q whose characteristic 
exponent is regularly varying in a special manner. 

Theorem 2.6. Choose and fix r G (0,2) and q G R. Then, there exists 
an isotropic Levy process X := {X t }t>o such that 

(2.55) < inf — - - < sup —. - ' „, < oo. 

^:||eir(l0g||CII) 9 -^R d: IICI| r (l0g||ell) 9 

U\\> e U\\>e 

Proof. Let B := {B(t)}t>o denote a d-dimensional Brownian motion, 
independent from the subordinator r := {rt}t>o of Lemma 2.5. Define 

(2.56) X t := B( Tt ) for all t > 0. 

It is well known — as well as easy to check — that the process X := {X t }t>o 
is a Levy process with characteristic exponent 

(2.57) E \e^ Xt ] = exp ( -t& (^-) for alH > and £ G R d . 



That is, *(£) = $(||£|| 2 /2); Lemma 2.5 [applied with p := r/2] completes 
the remainder of the proof. □ 



CHAPTER 3 



Positive-Definite Functions, Fourier Analysis, and 
Probabilistic Potential Theory 

A large part of this paper relies heavily on our ensuing analysis of 
positive-definite functions and their many connections to harmonic anal- 
ysis. In this chapter we develop the requisite theory and prove Theorem 
1.5. We also give intrinsically-probabilistic interpretations to the two cen- 
tral potential-theoretic hypotheses of this paper; namely Conditions 1.1 and 
1.10. Even though, most of the results derived in this chapter will be use 
later in this paper, they might be of independent interest. 

3.1. Fourier Analysis 

Let us begin by recalling some basic facts from harmonic analysis. 
First, let us recall our definition of Fourier transform. 

(3.1) g(0 := [ e ix< g{x) dx for all £ € R d and g G L l (R d ). 

This particular normalization is standard in probability theory and leads 
to the following form of the Parseval identity, which plays a big role in the 
ensuing theory: 

(3.2) / g( x ) h{x ) dx = J— f g(0W)^ forallg,heL 2 (R d ). 

The preceding is another way to say that the Fourier transform is an "isom- 
etry" from L 2 (R d ) onto itself. 1 

Recall that a function g : H d — > R is tempered if it is Borel-measurable 
and there exists k > such that 

(3 ' 3) i u £uTw <0 °' 



-^Strictly speaking, this is not true, as is evidenced by the multiplicative factor of (2n)~ d ^ 2 
in the right-hand side of equation (3.2). And that is why the word isometry appears in 
quotations: In order for the Fourier transform to be a proper isometry from L 2 (R d ) onto 
itself, we need to use a different normalization than the one used here. We have not done 
that as it would be nonstandard for probability theory. 
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Also, recall that g has at most polynomial growth [at least polynomial 
growth, resp.] if the preceding holds for some k > [k < 0, resp.]. Finally, we 
say that g G S if g and all of its derivatives have at least polynomial growth. 
The collection S is the usual space of rapidly- decreasing test functions on 

rA 

The following is an elementary, but important, variant of the Parseval 
identity. 

Lemma 3.1. Parseval' s identity (3.2) is valid when g G S and h : Fl d — > 
R is continuous and tempered. 

We leave the proof to the interested reader. 

3.2. Positive-Definite Functions 

Recall that a function g : H d — > R + is positive definite if g is tempered 
and (4>,g * <fi)L 2 (R, d ) > f° r an rapidly-decreasing test functions 4>. That is, 

(3.4) / (j){x)(j){y)g{x — y) dx dy > for all <f> £ S. 

JR d xIl d 

The following central result of L. Schwartz characterizes positive-definite 
functions. 

Proposition 3.2 (Schwartz [GV77, Theorem 3, p. 157]). If g : K d -»• R 

is positive definite, then there exists a tempered measure F on H d such that 
9 = t. 

The preceding has an elementary converse as well: Recall that a Borel 
measure T on H d is of positive type if F > in the sense of distributions; 
that is, r(0) > for all nonnegative (j) G S. Then it is straightforward to 
check that if F is a tempered measure of positive type on H d , then F is 
positive definite. 

Schwartz's theorem is a generalization of the following theorem of Her- 
glotz [d = 1} and Bochner [d > 2] : 

Proposition 3.3 (Herglotz [Herll], Bochner [Boc33]). If g : K d -)• R 

is continuous and positive definite, then there exists a finite Borel measure 
F on K d such that g = f . 

Recall also that a continuous function g : H d — > R is positive definite if 
and only if is is positive definite in the sense of Herglotz, Bochner, Polya, 
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etc.; that is if 

N 

(3.5) a i®j 9( x i ~ x j) > f° r au Oi,... ,ajv G C and x±, . . . ,xn G H d . 

The proof uses elementary integration theory, and we merely recall the [easy] 
steps: Define the finite complex measure 

N 

(3.6) :=J2 a i S xi, 

i=i 

and choose a sequence </>i, 02 , . . . : R d — > R of probability density functions, 
each in S, such that 4> n — > 1 pointwise as n — >■ oo. Then, one can make 
precise the following approximation: As n — > oo, 

N 

^2 a iOj9(xi - xj) = (g*fi)dfi 

(3.7) <J=i 

(g*(() n * fx)d((f) n * fx). 



And the final quantity is nonnegative because 5 is positive definite and 
Re(4> n * /i),lm(0 n * //) G 5 for all n > 1. 

3.3. A Preliminary Maximum Principle 

Now that we have recalled the basic definitions and properties of positive- 
definite functions, we can begin our proof of our maximum principle [Theo- 
rem 1.5]. But first let us prove the following technical result. 

Lemma 3.4. If (j) £ S, then there exists a version of f * <fi that is in 
C (R d ). Consequently, Rp(f * </>) G C (R d ) for every p > 0. 

Proof. Because / is tempered, the following defines a uniformly con- 
tinuous function on Tl d : 

(3.8) h(x) = J e~ ix <f(0H0 d£ for all x G R d . 

In fact, h G Co(R d ) because of the Riemann-Lebesgue lemma. Furthermore, 
if ip G S, then 



(3.9) 



/ ^(x)h{ X )dx = -± Yd I mKoko^ 

jR d i 2 ^) JR d 

= / i>{x)(f * 4>){x)dx. 
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The first line is justified by the Fubini theorem, and the second by the 
Parseval identity. It follows from density and the Lebesgue differentiation 
theorem that h = f * 4> almost everywhere. This proves the first assertion 
of the lemma. In addition, 



= {Rp(f *<!>)) (x). 

Since h G C (R d ), it follows from (2.10) that Rp(f*</>) = Rph G (7 (R d ). □ 



The following contains a portion of the said maximum principle of The- 
orem 1.5. It also provides some of the requisite technical estimates that are 
needed for the remainder of the proof of Theorem 1.5. 

Proposition 3.5. For all (3 > 0, 



(3.11) T(/3) = sup (Rpf)(x) = ess 8up(Rpf)(x) = limsup( J R /3 /)(x), 

a;GR d x£R d x^>0 

where T(0) := lim^ T(/3). 



Proof. First, we prove the proposition in the case that /3 > 0. 
In accord with Lemma 3.4, if (f) G S, then Rg(f * (ft) is continuous. 
Therefore, the Plancherel theorem applies pointwise: For all x G R rf , 



[Without the asserted continuity, we could only deduce this for almost every 
x G R d .] In particular, for all probability densities <j) G S, 

(3.13) sup (Rpif *</>)) (x) <T(/3). 

[T was defined in (1.15).] If {(p n }^ =l is an approximate identity consisting 
solely of probability densities in S, then 

(3.14) Rgf < liminf Raff * 4> n ) [pointwise], 

n— >oo 

by Fatou's lemma. Consequently, (3.13) implies that 



(3.10) 




(3.12) 




(3.15) 



sup (R fi f)(x)<T(P). 
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In order to prove the reverse bound, define the Gaussian mollifiers {4> n } 



oo 

n=l- 



I 112 

\z\\ n 



(3.16) 4> n {z):= (|-) d/2 exp(- 

And observe that 

If f7f) e -|l«H 2 /(2n) 

(317) ri( 

= (l + o(l))T(/3) as n ^ oo, 
thanks to the monotone convergence theorem. On the other hand, 

{Rp(f*<f>n))(0)= [ rf,(y)(f*<f> n )(y)dy 



(3-18) = / {Rf,f)(y)Mv) 



dy 



< ess sup(Rpf)(x), 

x£~R d 



since J Rd (f) n (y) dy = 1. This and (3.13) together prove that T(/3) is the 
maximum of the /3-potential of /; and the maximum ^-potential is finite 
if and only if T(/3) is. We can choose the </> n 's so that in addition to the 
preceding regularity criteria, every 4> n is supported in the ball of radius 1/n 
about the origin. In that way we obtain 

{Rpif*^)) (0) = (W)*^)(0) 

(3 ' 19) < sup {Rpf)(x). 

\\x\\<l/n 

And this proves (3.11) for every /3 G (0 , oo). 

The case that (3 := has to be handled separately because is not a 
finite measure. Therefore, we next derive (3.11) in the case that (3 = 0. 

First of all, j3 i->- T(/3) is nonincreasing. Therefore, T(0) := lim^o T(/3) 
exists as a nondecreasing limit. Because Rpf < Rof [pointwise] for all 
(3 > 0, we can deduce that 

(3.20) T(0) < ess sup {R f )(x), 

x£R d 

and 

(3.21) T(0) < limsup(i? /)(x). 
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For the reverse bound, recall (3.13), and let ft \. 0. Since T(/3) — > T(0), 
we find that for every probability density (f> G S and x G R ', 



(3.22) 

But the left-hand side is 



ljm(iW*<M) (*)<T(0). 



(3.23) 



limE 



(f * ^(Xs + x)e-? s ds 



thanks to the monotone convergence theorem. Another application of Fa- 
tou's lemma shows that (Rof)(x) < T(0) for all x G R d . This establishes 
(3.11), and hence the proposition, in the case that ft = 0. □ 



3.4. Proof of the Maximum Principle 

The main goal of this subsection is to establish Theorem 1.5. This sub- 
section also contains a harmonic-analytic estimate that might be of indepen- 
dent interest. We will use that harmonic- analytic to demonstrate Theorem 
1.5, as well as the subsequent Theorems 1.8 and 1.11. 

In order to motivate our estimate, let us first consider the important 
special case that the correlation function / is of Riesz type. That is, 

(3.24) f(z) := \\z\\-( d -V for z G R d , 

where _1 := oo. Clearly, / is locally integrable when b G (0 , d), a condition 
which we assume, and in fact its Fourier transform is 

See Mattila [Mat95, eq. (12.10), p. 161], for example. Then, it is well known 
that 

H(dx)n(dy) 



J J f(x-y) f i(dx) f i(dy) = jj 



f - y\\ 

|2 



(3.26) = *w*r m r jacol 

r((d-6)/2) 7 Rd neii 6 ? 

1 * mwnodt. 



(2vr) d 7 Rd 

It is easy to guess this famous identity from an informal application of the 
Fubini theorem. However, a rigorous derivation of (3.26) requires a good 
deal of effort; see Mattila [Mat95, Lemma 12.12, p. 162], for instance. 
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In the language of potential theory, the preceding asserts that the Riesz- 
type "energy" of the form ff \\x — y\\~ d+b n(dx) n(dy) is equal to a constant 
multiple of the Polya-Szego-type energy J Rd |A(£)| 2 |l£ll~ fe d£. The following 
proposition shows that there is a very general lower bound that is valid for 
every correlation function /. 

PROPOSITION 3.6. For all Borel probability measures \i on H d , 

(3-27) J J f{x - y) fi(dx) M (dj/) > J^y f^ d IA(0l7(0 d£ 

Proof. We will assume, without incurring any loss in generality, that 

(3.28) J J fix - y) nidx) n(dy) < oo; 

for there is nothing to prove otherwise. We also observe that 

(3.29) ff fix - y) M (ds) ti(dy) = J (f * dp. 

This is essentially the famous "reciprocity theorem" of classical potential 
theory. 

Because / * ll € L 1 ^), Lusin's theorem implies that for all e > there 
exists a compact set A e in R rf such that: 

(i) n(A°) < e; and 

(ii) / * fx is continuous on A e . 

Let 

(3.30) Me(*) := fi{»nA e ) 



denote the restriction of /j, to A e , and recall the Gaussian densities {<f\ 
from (3.16). Because 

(3.31) lim (/ * 4> n/2 * n) = f * fx, uniformly on A e , 

it follows from (3.29) that 

ff fix- y) nidx) n(dy) > f if * fj) d^ 

(/ * n) d/i e 



oo 
n=l 



(3.32) 

lim / (/ * 4> n/2 * n) dfx e 
"~> J 



n— >oo 



> limsup / (/ * <p n/2 * fJ, e ) d/J, e . 

n—s-co J 
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(3.34) 



Since <f> n / 2 = 4> n * 4>n, 

(3.33) / / f(x- y) n{dx) n(dy) > limsup / (/ * /v e )(z)// n , e (x) dx - 

where fi n;e (x) := {4> n * (i e )(x). Since ^ n)£ G S and / * ^ ni£ is C°° and 
tempered, Lemma 3.1 implies that 

/ (f*Hn,e)(x)^n,e( x ) dx =-pr-Td[ /(OIAn, e (0| 2 

=(2^L eHI£llv " i ''<« )|2/ «» d «- 

This, (3.32), and the monotone convergence theorem together imply that 
(3-35) J J f(x - y) ^(dx) M (dy) > -^p |A e (0l7(0 d£- 

Since ^ e converges weakly to /i as e — > 0, we know that /t e converges to /t 
pointwise. The proposition follows from this and Fatou's lemma. □ 

An elementary computation [Mat95, Lemma 12.11, p. 161] and Propo- 
sition 3.6 above together yield the following corollary. 

Corollary 3.7. // / is lower semicontinuous, then for all Borel prob- 
ability measures ji on H d , 

(3-36) JI f(x - y) M (dtf) /i(dy) = J^yi J IA(0l7(0 d£. 

Even though we do not use the above result in the sequel, we have chosen 
to include it as it is an interesting fact about energy forms that are based on 
lower semicontinuous positive-definite functions. A weaker version of this 
result has been proved by Khoshnevisan and Xiao [KX09, Theorem 5.2]. 

We are now ready to prove Theorem 1.5. 

Proof of Theorem 1.5. Note that for all t > 0, 

(p t f)(o) = (p t *p t f) (o) 

(3 >37 ) = (pt *Pt*f) (0) 

= H f(x-y)pt(x)p t (y)dxdy. 

This requires only the Tonelli theorem. Now Proposition 3.6 can be applied 
to show us that 

(3.38) (A/)(0) >J^yij nd e- 2tRe * ( «/(6 d*. 
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Multiply both sides by exp(— fit) and integrate [dt] to find that 

(3.39) ( ^ /)(0) - wr L i+iBm de 

= T(/3). 

This and Proposition 3.5 together imply (1.17). Consequently, (1.1) holds 
if and only if T(/3) < oo for all /3 > 0. 

Next we prove that Rpf G Co(R d ) whenever T(/3) < oo and / is lower 
semicontinuous. 

When / is lower semicontinuous we can find compactly-supported con- 
tinuous functions f n that converge upward to / as n — > oo. Recall from 
(2.10) on page 15 that Rp maps Co(R d ) to Co(R rf ). Consequently, Rpf n G 
Co(R d ), and from this we may conclude that Rpf is lower semicontinuous. 

Next, let us define 

(3.40, Ml):= _^_m. de 

If T(/3) < oo, then irp G Co(R d ). Moreover, a few successive applications of 
Fubini's theorem tell us that for all <fi € S, 

(3.41) / Tip{x)4>{x)dx= f {Rpf){x)4>{x)dx. 

JR d JR d 

Thus, iTp = Rpf a.e. 

It remains to prove that if T(l) is finite, then so is T(/3) for every j3 > 0. 

We have shown that if Y(l) is finite, then R\f = tt\ almost everywhere, 
7Ti G Co(R d ); also, R\f is bounded [Proposition 3.5]. If h\ = h 2 almost 
everwhere and hi,h 2 : R d — >■ R+ are measurable, then 

(Rph 1 )(x) = fp{y)h 1 {x -y)dy 
Jn d 

(3 - 42) = / fp(y)h 2 (x-y)dy 

Jn d 

= (Rph 2 )(x), 

for all x G K d . Therefore, (2.10) on page 15 and the fact that 7Ti G C (R d ) 
together imply that Rpir\ G C'c^R ') for all /3 > 0. In particular, RpR\f G 
Co(R d ) — whence RpRif is bounded — for every /3 > 0. And by the resolvent 
equation [see Blumenthal and Getoor [BG68, (8.10), p. 41]], 

(3.43) (Rpf)(x) = (R l f)(x) + (l-(3)(RpR l f)(x) for all x G R d . 
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Thus, it follows from the boundedness of R±f and RpR\f [for all f3 > 0] 
that Rpf is bounded for every (3 > 0. Consequently, Proposition 3.5 implies 
the finiteness of T(/3) for every (3 > 0. □ 



3.5. Probabilistic Potential Theory 

Define, for all measurable functions 4> : H d — > R + , a process t i-» L t ((f)) 
as follows: 

(3.44) L t ((j)):= f 0{X s )ds for alH € [0 , oo] . 

Jo 

The random field (t,<f>) >->■ Lt(</>) defined above is often called the occupa- 
tion field of the Levy process X; see, for example, Dynkin [Dyn84]. It is 
well defined, but might well be infinite even in simple cases. The following 
example highlights this, and also paves the way to the ensuing discussion 
which yields a probabilistic interpretation of Conditions 1.1 and 1.10. 

Example 3.8 (After Girsanov, 1962 [McK05, §3.10, pp. 78-81]). Con- 
sider d = 1, and let X denote one-dimensional Brownian motion, normalized 
so that 

(3.45) Ee^' Xt = e"* f2/4 for alH > and (El 

In this way, X is normalized to be standard linear Brownian motion; that 
is, 

(3.46) Ee i? '^ = e"* f2/2 for all t > and £ G R. 

Let <j){x) := |x|~ Q for all x G R, where a G (0 , 1) is fixed. By the occupation 
density formula [see Corollary (1.6) of Revuz and Yor [RY91, p. 209]], 

/oo ex 
-oo 

where £ denotes the process of local times associated to X. According to 
Trotter's theorem [RY91, Theorem (1.7), p. 209] and the occupation density 
formula, 

(3.48) P {£' G C (R) for all t > 0} = 1. 

It is a well-known consequence of the Blumenthal zero-one law and Brownian 
scaling that 

(3.49) P {4° > for all t > 0} = 1. 
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Consequently [if we ignore the null sets in the usual way], L t {4>) < oo for 
some — hence all — t > if and only if a < 1. At the same time, we note that 

rt ds 



EL t {^>) = E 

(3.50) 



/ 

J o 



Jo W 2 V^J-oo 

This requires only Brownian scaling and the Tonelli theorem. Thus, EL t ((f)) 
is finite for all t > 0, if and only if a < 1. In rough terms, we have 
shown that Lt{4>) < oo if and only if ELt{<p) < oo. As we shall see, a 
suitable interpretation of this property can be generalized; see Theorem 
3.13 below. □ 

It is convenient to use some notation from Markov-process theory: Recall 
from Markov-process theory that P z denotes the law of the underlying Levy 
process started at z 6 R rf [so that P = Po], and E z denotes the corresponding 
expectation operator. Since the underlying Levy process X is Levy, V z can 
be interpretted as the law of X, + z. Thus, 

poo 

(3.51) (R a <f>)(z)=E g e- as <P(X s )ds. 

Jo 

Before we state and prove the main result of this section, let us first 
establish some technical facts which will be needed in the proof of the main 
result. 

Lemma 3.9. For all t, a > and measurable functions 4> '■ R d — > R-+, 

(3.52) sup E x (sup [e" as L s (0)] ) < sup (R a 4*){z) < Xa (t) sup E x L t {4>), 

xe~R d \s>0 J zG~R d x£R d 



where 

(3.53) X a(t) := 



e at _ 1 



Remark 3.10. Let us point out the following elementary bound for the 
right-most term in the preceding display: For all t, a > 0, 

(3.54) e~ at sup E x L t {4>) < sup E x (sup [e" as L s (0)] ) . 

xeR d xeR d \s>o / 

This shows that the quantities on the two extreme ends of (3.52) are one 
and the same, up to a multiplicative constant that depends on a and t in 
an explicit way. □ 
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Proof of Lemma 3.9. Because (j) is nonnegative, we have the sure in- 
equality, 

(3.55) e~ as L s ((P)< e~ ar (p(X r )dr, valid for all s, a > 0. 

Jo 

We take suprema over s > and then apply expectations [dP x ] to deduce 
the first inequality in (3.52). For the second bound, let us note that for all 
a, t > and z G R d , 

00 r(n+l)t 

(R a( f>)(z) = J2 E z e- QS 0(X)d S 

n J nt 



n=0 



(3.56) 



n=0 

This implies the second inequality, because 

ft rt 



(3.57) 



E z / <P{X s+nt ) ds = E z E Xnt [ <f>(X a )da 
Jo Jo 



< sup [L t 
xeR d 



in accord with the Markov property. □ 
Lemma 3.11. For all t > and measurable cf> : R d — > R + , 

(3.58) sup (\L t {4>)\ 2 ) < 2 ( sup E z [L t (4>)] \ . 

Proof. We can write 

(3.59) E x (|L t (0)| 2 ) = 2 jTdu jf dv E, (0(X U )^(X)) • 
Since 



(3.60) 



E^ ^ 0(X) dw X s ; s< u ^j =E Xu j^ <j>{X 8 ) ds 

= E Xu [L t - U {4>)} P x -a.s., 



it follows that 



(3.61) 



E* (|^(0)| 2 ) = 2E X j%{X u )E Xu [L t _ u (4>)} du 

< 2E, / 0(X U ) • sup E, [L t _ u (0)] du. 

JO zgRd 



The lemma follows since Lt_ M (<^>) < L t ((j)). □ 
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The following constitutes the third, and final, technical lemma of this 
section. 



Lemma 3.12. For all a, t > 0, 

(3.62) s W L ^%^r}4 

where Xa(t) is defined in Lemma 3.9. 



PROOF. Recall the Paley-Zygmund inequality [PZ32]: If Z is a nonneg- 
ative random variable in L 2 (P) with EZ > 0, then 

(3.63) p{z>1 E ,}>0-. 

We can apply this with Z := L t (4>) — where 4> : H d — > R + is bounded away 
from zero and infinity — to see that for all z € R rf , 

(3.64) sup vjl.no; > ii;. /,ioil • ( ^ L ' (o)) ' 



By selecting z appropriately, we can ensure that 
(3.65) E z L t {4>) > (1 - e) • sup E x L t (<f>), 

x£R d 

where e G (0 , 1) is arbitrary but fixed. Let e | and appeal to the continuity 
properties of probability measures to deduce that 



(3.66) 



D J r a\ n 1 n r^L ( su P*eR d V z L t {4>)) 2 
sup P x < L t (4>) > - sup E z L t {(p) } > 



xeR" [ 2 Zen* J 4sup :rgRd E ;c (|Li((/))| 2 ) 

1 

> -; 



see Lemma 3.11 for the last inequality. A monotone-class argument shows 
that the preceding holds true for all bounded and measurable functions 
4> ^ 0. We apply it with 4>n := min(/ , N) in place of 0, where N > 1 is 
fixed. In this way we obtain 

(3.67) sup P x \ UU) > I sup E.L^n) \ > 1. 

x£~R, d ! Z z€R d 1 5 
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This and Lemma 3.9 together tell us that 

sup p x U t (f) > iMi^M ) 

( 3 - 68 ) >su P p4m/)>su P ^%^M) 

xe~R d { zeR d z Xa\t) J 

1 

> -. 

~ 8 

As N t oo, (R a (p N )(0) t (R a f)(0), and the lemma follows. □ 

The next result yields a probabilistic characterization of Condition 1.1 
which we recall for the reader's convenience. 

{R a f)(0) < oo for all a > 0. 

Theorem 3.13. Under Condition 1.1, 

(3.69) P z {L t (f) < oo for all t > 0} = 1 for all z G R rf . 
Moreover, in this case, t >->■ L t (/) grows subexponentially. That is, 

(3.70) P z (limsup log -W) < j = 1 for all ze H d . 

{ t^QO t J 

On the other hand, if Condition 1.1 fails to hold, then 

(3.71) P 2 {L t (f) < oo for some t > 0} = /or some z £ R rf . 

Remark 3.14. Consider the stochastic heat equation where F is space- 
time white noise. Formally speaking, this means that / := Sq is our corre- 
lation "function." In this case, one can [again formally] interpret 

(3.72) L t (f) = L t (S ) = /\(* s )ds 

as the local time of the replica process X at zero. And if we interpret 
Theorem 3.13 loosely as well, then Theorem 1.2 suggests that (1.1) has a 
mild solution if and only if X has local times. This interpretation is correct, 
as well as easy to check, and leads to deeper connections between SPDEs 
driven by space-time white noise on one hand and local-time theory on the 
other hand [EFK09, FKN09]. In the case of the parabolic Anderson model 
[that is, (1.1) with cr(u) = const - u and 6 = 0], Bertini and Cancrini [BC95] 
and Hu and Nualart [HN09] discuss other closely-related connections to 
local times. In the case that £ is a discrete variable [for example, because C 
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is the generator of a Levy process on Z d ; i.e., the generator of a continuous- 
time random walk], similar connections were found earlier; see Carmona and 
Molchanov [CM94], for instance. □ 

Proof of Theorem 3.13. If (R a f)(0) < oo for some a > 0, then 
(Rpf)(0) < oo for all /3 > by Theorem 1.5. It follows from the first 
inequality of Lemma 3.9, and Theorem 1.5, that for all (3 > and z G R rf , 



(3.73) E 



sup 

t>0 



< (Rpf)(0) < oo. 



(e"*M/)) 
This implies (3.69); it also implies that 

(3.74) limsup e~ tf5 L t (f) < oo almost surely [P 2 ]. 

t— >oo L 

This implies (3.70) because (5 > and z G R d are arbitrary. 

In order to finish the proof, let us consider the remaining case that 
CRa/)(0) = oo for all a > 0. 

According to Lemma 3.12, 

(3.75) sup P^ {L t (f) = oo} > - for all t > 0. 

In particular, there exists z G R d such that 

(3.76) P 2 {L t (f) = oo for some t > 0} > -. 

Because t >->■ L t (/) is nondecreasing, the Blumenthal zero-one law applies 
and implies that 

(3.77) P z {L t (f) = oo for some t > 0} = 1; 

this implies the remaining portion of the theorem. □ 

We now have the following consequence of Theorem 3.13. It is particu- 
larly useful because its hypothesis are verified by all the examples that we 
have mentioned in the Introduction. 

Corollary 3.15. Suppose f is bounded uniformly on the complement 
of every open neighborhood of the origin. Then, Condition 1.1 is equivalent 
to the following: P{Lt(f) < oo for some t > 0} = 1. 

Proof. According to Theorem 3.13, if Condition 1.1 holds then 

(3.78) P z {L t (f) < oo for some t > 0} = 1 for all z G K d . 
Set z := to obtain half of the corollary. 
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Conversely, suppose Condition 1.1 fails. According to Theorem 3.13, 
there exists a point z G R d such that 

(3.79) P z {£*(/) < oo for some t > 0} = 0. 

We need to prove that z = 0. This holds because if z were not equal to the 
origin, then 

(3.80) P z {L t (f) < oo for all t € [0 , t)} = 1, 

where r denotes the first hitting time of the open ball of radius \\z\\/2 around 
0. Indeed, 

(3.81) sup L t (f) < t ■ sup f(u) < oo, 

0<t<r \\u\\>\\z\\/2 

P z -almost surely. Since the paths of X are right-continuous, P z {r > 0} = 1, 
and hence (3.79) is contradicted. □ 

Condition 1.10 [p. 8] also has a probabilistic interpretation that is given 
by following proposition. 

Proposition 3.16. If (R a f)(0) < oo for some, hence all, a > 0, then: 

(3.82) (Rof)(0) < oo => Loo(/) < oo a.s.; 
and 

(3.83) ( J R /)(0) = oo Loo(/) = oo a.s. 

Proof. If (Rof)(0) < oo, then because (^o/)(0) = EL^/) < oo, it 
follows that L OQ (f) is finite a.s. 

If, on the other hand, (_Rq/)(0) = oo, then because 

poo 

(3.84) E/ e- QS /(X)d S = (i? a /)(0), 

JO 

the Paley-Zygmund inequality (3.63) implies that 

(3.85) p{ re- as f(X s )ds > l(R a f)(0)) > |(^/)(0)| 



4E(|/ °°e— f(X s )ds\ 2 )' 



It follows from this and Lemma 3.11 — see also the proof of Lemma 3.12- 
that 

p{ rnx.)** > 1(rmo)\ > - ( ^ 0) nM 

(3. 86 ) Uo 2 J 8sup a;gR d( J R a /)(x) 
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owing to Theorem 1.5 [p. 5]. Let a I to find that 

V:=pIJ^ f{X s )ds = ^ 



(3.87) 

> -. 



But according to Theorem 3.13, J Q T }'{X S ) ds < oo, for all T > a.s., because 
(R a f)(0) < oo for some [hence all] a > 0. This implies that 

(3.88) V = P^lim j f(X s )ds = oo\. 

That is: (i) V is the probability of a tail event; and (ii) V is strictly positive, 
in fact V > 1/8. By the Hewitt-Savage zero-one law [HS55], V = 1. □ 



3.6. A Final Observation 

Let us conclude this chapter with an observation that will be used later 
on in Theorem 4.9 [p. 51] in order to produce a stochastic PDE whose 
random-field solution exists but is discontinuous densely. 

Let X := {X t }t>o denote a Levy process on R d with characteristic 
exponent 'I'. Recall that X has a one-potential density v if v is a probability 
density on H d that satisfies the following for all Borel-measurable functions 
6 : R d -»• R+: 



f'OO P 

j e~ s (p(X s )ds = / cp(x)v(x)dx. 
JO JR, d 



(3.89) E 

Because 4> > 0, the preceding expectation commutes with the ds-integral. 
Recall that m s denotes the law of X s , and restrict attention to only non- 
negative (f> € S. In that case, 

/'OO f'OO / f \ 

J e~ s E(P(X s ) ds = J e" s (j <f>dm 3 \ ds 

(3.90) -&?r'"_*f» i(M '~™ 

1 f *«> ae 



(2»>j y RJ i + *({) 

We compare this to the right-hand side of (3.89), and then apply Plancherel's 
theorem to the latter, to deduce the following well-known formula: 

(3.91) v(0 = for all ^ G R d . 
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If we consider only the case that X is symmetric, then v is rendered 
nonnegative, since ^ is nonnegative in this case. This observation and the 
Bochner-Schwartz theorem [Theorem 3.2] together imply that v is a correla- 
tion function. Because products — and hence integer powers — of correlation 
functions are themselves correlation functions, Theorem 2.6 yields the fol- 
lowing byproduct. 

Theorem 3.17. Choose and fix a > and b G R. Then, there exists a 
correlation function v on H d such that 

(3-92) v(0 ^ for e € R d wtth U\\ > e. 



CHAPTER 4 



The Linear Equation 



Before we study the fully nonlinear equation (1.1), we analyse the far 
simpler linearized form of the same equation [(7=1, 6 = 0], and show that 
it has many interesting features of its own. Because the solutions, if any, to 
the said linear equations can only be Gaussian random fields, we are able 
to use the theory of Gaussian processes in order to produce some definitive 
existence and regularity results. Our results should be compared with our 
earlier joint effort with Eulalia Nualart [FKN09], in which F was space- 
time white noise. Our earlier effort was, in turn, motivated strongly by the 
earlier works of Dalang and Frangos [DF98], Dalang [Dal99], and Peszat 
and Zabcyzk [PZOO]. 



subject to uq being the initial function, which as mentioned in the Intro- 
duction is assumed to be a nonrandom bounded and measurable function 
uq : H d — > R. One can follow through the theory of Walsh, and define the 
weak solution to (4.1) as the Gaussian random field u := {«<(<?!>) }t>o,</>eS, 
where [we recall] S denotes the collection of all rapidly-decreasing test func- 
tions on R d , and 



The double integral is a Wiener integral and {Pf}t>o is the semigroup as- 
sociated to the dual process X t * := — X t [t > 0]. 

There are two main questions that one needs to answer before one pro- 
ceeds further: 

(a) Is u well defined? 

(b) What is the largest family of 0's for which ut(4>) is well defined? 



4.1. Existence and uniqueness 



(4.1) 



The linearized form of (1.1) is the stochastic PDE 
) ^-u t (x) = (Cu t )(x) + F t (x), 



(4.2) 
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An affirmative answer to the first question would imply existence of solutions 
in the general sense of Walsh [Wal86]. Since the analysis of the nonrandom 
quantity J Rd uq(x)(P^ 4>){x) dx is standard, we can reduce our problem to 
the special case that uq = 0. In that case, these question are addressed by 
the following estimate. Here and throughout, we define 

(4 - 3) £ ^--=w^L^T^m d ^ 

for all Schwartz distributions v whose Fourier transform is a function. 

Lemma 4.1. The weak solution u to (4.1) with uo = exists as a well- 
defined Gaussian random field parametrized by t > and <p € S. Moreover, 
for all t, A > and <f> 6 S, 

(4.4) a(t)S x (<t>) < E (|n t (0)| 2 ) < b{t)£ x {<t>), 
where a{t) := (1 - e~ 2t l x ) and b{t) := e 2 */ A . 

In fact, Lemma 4.1 holds under far greater generality than the one pre- 
sented here. For instance, it holds even when transition functions do not 
necessarily exist, and when the correlation function is a general correlation 
measure. The less general formulation above suffices for our needs. 

Proof. If <f> G S then pt<p G S for all t > 0. Since the Fourier transform 
is an isometry on S, this proves that 

(4.5) Pl<\> = p t * (f> G S for every t > 0. 

Therefore, in accord with (1.3), the second moment E(\ut(4>)\ 2 ) is equal to 



E 



f I (P t *- s 4>)(y)F(dsdy) 
Jo Jn d 



2^ 



1 /•* 



• worm 



(2vr)^y y Rd i 

= -J— f\, f <\P p-2»Re*(0 . 
(2vr) d 

And the lemma follows from the preceding and Lemma 3.5 of [FKN09]. □ 

There are standard ways to extend the domain of Gaussian random 
fields. In our case, we proceed as follows: Consider the pseudo-distances 
{pt}t>o defined by 

1 /2 

(4.7) Pt{<f>,ip) ■= {e(m</>)-^(vO| 2 )} for </>,Ve<s. 
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Because L 2 (P)-limits of Gaussian random fields are themselves Gaussian 
random fields, we deduce the following: Suppose v is a Schwartz distribution 
such that lim^oo p t (v , v * 4> n ) = for all t > 0, where {<fi n }^ =1 is the 
sequence of Gaussian densities, as defined in (3.16) [p. 29]. Then u t (v) is 
well defined in L 2 (P), and the totality {ut(v)} of all such random variables 
forms a Gaussian random field. 

We follow [FKN09] and say that (4.1) has a random-field solution if we 
can obtain Ut(5 x ) in this way for all x G R rf . 

Consider the space %$ of all <p G S such that £\{4>) < oo for all t > 0. 
Evidently, %q can be metrized, using the distance 

(4.8) S(<f>,fp) := v^iOA-VO forces. 

Define %i to be the completion of Ho in the distance 5. 

Lemma 4.2. Condition 1.1 holds iffS x G %i /or some, hence all, x G R d . 



Proof. First of all, we recall (1.15) and check that 

(49) £x{Sx) = W*T d L A- 1 + 2Re^(0 

= ^T(l/A). 

In particular, the value of £\(5 X ) does not depend on x G R d . And Theorem 
1.5 implies that £\(5 X ) is finite for some A > if and only if it is finite for 
all A > 0. 

Let us first suppose that £\(S X ) is finite. We can note that S x * 4> n = 
cj) n {» — x) G S, where {(p n }'^L 1 was defined in (3.16). Therefore, for all 
n, m > 1, 

£\(&x *4>n-$x* 4>m) 



(4.10) i r i 



1 1 2 I 1 1 



\ _ g llsll | 2n 2m i 



2 



2(2vr) d 7 Rd A- 1 + 2Re*(0 

Since £"a(5o:) is finite, the dominated convergence theorem tells us that the 
sequence {S x *4> n }^ =1 is Cauchy in %q. A calculation similar to the preceding 
shows that the quantity £\(5 X * (j) n — S x ) converges to zero as n — > oo. And 
therefore, 5 X G Hi. 

Conversely, if S x G Hi, then £\{8 X * 4> n — S x * (j) m ) — > as n, m — > oo. We 
can extract an unbounded subsequence n± < n<i < ■ ■ ■ of positive integers 
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such that 

(4.11) £ X {5 X * </> nj - S x * cp nj+1 ) < 2-' for all j > 1. 

It follows from (4.10) that if k~ l < {nj 1 - nj^], then 

£x(8 x - S x * cf> k ) = - JL, f — 1 - e-^ll 2 /^) 

"'27r) d 7 Rd A- 1 + 2Re*(C) 



(412) <S x (6 x *<l> ni -S x * 



2(27 



2 



< 2^'. 

Let — > oo and then j — > oo, in this order, to deduce from the preceding 
discussion that 

(4.13) lim S X (6 X - 6 X * <f> k ) = 0. 



Because v \-> y / £\(v) satisfies the triangle inequality, it follows from (4.11) 
that for all k > 2, 



y£>jjx) < ythJJx - S x * (f> k ) + y/£\{5 x * <Pk) 



< yf£\(5 x -$x* <t>k) + V £>JJx *4>k-$x* <t>k-i) 

+ \/£\{$x * 4>k-i) 

(4.14) < V£Wx - $x * A) + 2-( fc - 1 )/ 2 + y/£ x (6 x * <^-i) 



fc-i 

< V£\(6x-6x*<l>k) + 2 ~ j/2 + y/£x(s x *</>!). 

3=1 

We have shown that the first quantity on the right-hand side converges to 
zero as k — > oo; and the second term remains bounded. Finally, the third 
quantity on the right-hand side of the preceding is finite since <j)\ G S. 
Therefore, it follows that if S x G Hi then £\(5 X ) < oo. This concludes our 
proof. □ 

For more general initial functions uq > 0, (4.1) has a random-field solu- 
tion if and only if Condition 1.1 holds and (P t uo)(x) < oo for all t > and 
x G R rf . Let us conclude this section with a lemma that provides simple 
conditions that ensure that (PtUo)(x) is finite for all t > and x G H d . 



Lemma 4.3. Suppose exp(-Re^) G L*(R d ) for all t > 0, and ii G 
//(R 6 *) /or some f3 G [l,oo]. T/ien, (P t u )(x) < oo /or all t > and 
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x G R d . Moreover, PtUo is uniformly bounded and continuous for every 
fixed t > 0. 

PROOF. Since Pt is a contraction on L°°(R d ), it suffices to consider only 
the case that 1 < (3 < oo. 

Choose and fix some t > 0. According to Young's inequality, 

\\ P tUo\\ L ^md) = \\pt*uo\\L°°(R d ) 
(4.15) 1 ' ( ' 

< l|Pt|lLP(R d ) - ll u ollL9(Rd), 
where p' 1 + c/ _1 = 1. On the other hand, for all p G (1 , oo), 

( 4 -!6) IHIlp(R<*) < IMIl°°(R<0> 

and this is finite, thanks to Proposition 2.3 on page 19. Therefore, it remains 
to prove continuity. 

First consider the case that ft < oo. In that case, we can bound the 
quantity 

(4.17) \(P t u )(x) - (P t u )(x')\ = \(p t *u )(x) - (p t *u )(x')\ , 

from above, by 

Pt(y) \uo{y - x) - u (y - x')\ dy 

VP 



L 



/R d 

(418) < ^ d Pt(y)\uo(y-x)-u (y-x')f dyj 

< ||Pt|lii (Rd) • \\M* ~ x ) ~ M* ~ x ')\\ L f3(nd) ■ 

It is a classical fact that u$ G Z/(R d ) implies that uq is continuous in 
L /3 (R d ). Therefore, PtUQ = p t *uo is continuous. 

Next let us consider the case that (3 = oo. In that case, we write 

\(p t *u )(x) - (p t *u )(x')\ 
(4.19) , 

< ll«0||L->(Rci) • \\Pt(» - X) - p t (» - X )\\ Ll(Rd) , 

which goes to zero because, once again, p t € L x (R d ) implies that pt is 
continuous in L 1 (R d ). □ 

4.2. Spatial regularity: Examples 

Define for all x, y G R ', 



(4 , 0) ^■^(^/ a /-~^-'' )) /(Od f ) 



1/2 
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Then, d defines a pseudo-distance on R d . Let Nd denote the metric entropy 
of [0,l] d . That is, for all e > 0, N d (e) denotes the minimum number of 
radius-e d-balls required to cover [0, l] d . We can combine Lemma 4.1 with 
theorems of Dudley (see, for example, Marcus and Rosen [MR06, Theo- 
rem 6.1.2, p. 245]) and Fernique [MR06, Theorem 6.2.2, p. 251], together 
with Belyaev's dichotomy [MR06, Theorem 5.3.10, p. 213], and deduce the 
following: 

Proposition 4.4. Suppose exp(-Re^) G L*(R d ) for all t > 0, and 
T(l) < oo so that (4.1) has a random-field solution u with uq = 0. Then 
the following are equivalent: 

(1) x I—?- ut(x) has a continuous modification for some t > 0; 

(2) x i y ut{x) has a continuous modification for all t > 0; 

(3) The following metric- entropy condition holds: 



Next we describe a large family of examples for (4.1) that have continu- 
ous random-field solutions. Throughout, we write "/i x g" in place of ll (h/g) 
is bounded, above and below uniformly, by finite positive constants." 

Theorem 4.5. Suppose f(x) = const/ and Re^(^) x ||£|| Q for 
some a £ [0,2] and /3 £ (0,<i). Then (4.1) has a random-field solution if 
and only if a + j3 > d. In this case, (4.21) holds. In fact, we have 

(4.22) d(x ,y) >z g(\\x — y\\) uniformly when \\x — y\\ < 1/e, 

where for all r G (0 , l/e) ; 



Remark 4.6. The condition that Re^(£) x ||£|| a implies that the upper 
and lower Blumenthal-Getoor indices of ^ match and are both equal to 
a; see Blumenthal and Getoor [BG61, Theorem 3.2] and Khoshnevisan 
and Xiao [KX09] for definitions and further details, including the various 
connections that exist between those indices and the fractal properties of 
the underlying Levy process X. □ 

Proof of Theorem 4.5. Recall (1.15). In order to prove the existence 
of random-field solutions, it suffices to show that T(l) < oo. We begin by 



(4.21) 




(4.23) 



g(r) := < r y / log(l/r) ifa + /3 = d + 2, 
r ifa + (3>d + 2. 
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writing 

T(l) = - J- / f{0 
(4.24) V ; (2vr) d J Rd 1 + 2Re*(£) 

:=/i+/ 2 , 



where 



1 f /«) d£, and 



(4.25) 



• (2vr) d 4||< 1 l + 2Re^(e) 

r - 1 / /(O j<- 

2 - (27T)" ^iim 1 + 2Ref (0 

It is clear from the hypothesis of the Lemma that I\ is always finite, because 
(3 < d. We now turn our attention to I2, and note that 

( 426 ) h x / TT^- 

J||?||>i IK II p 

Therefore, 22 < °o if and only if a + (3 > d. This concludes the first part of 
the result. 

For the second part we assume that 

(4.27) \\x-v\\ < 1/e. 

We write 

(4-28) | d(X)y )|2 := _L_ ( j 1 + J 2 + J 3 ), 



(2vr) c 



where 



W 1 7°f R ( ;- j ' )) /K)^. 

•%ll<i l + 2Re*(£) 

< 4 - 29 > W 1 7°f Re t7 )) /K)d ( , 

•%ll>VII*-2/ll l + 2KeV(^j 
Jl<lieil<l/lb-y|| l + 2KeV(4J 



L<||£||<l/||s-j, 
We can estimate each Jj separately. 
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Because 1 - cos 6 x 9 2 for 9 G (-1 , 1), 

Ji^\\x-yf- I ||el| 2 /(0d£ 

hs\\<i 

-> . /" lltl|2-/9, 

£11 

1 dr 



(4.30) >c ■ / iieir^ 



|x - 2/ 1 1 2 



Because (3 < d, the integral term in the above display is finite, and hence 

(4.31) J\ x ||x - y|| 2 . 

We estimate J2 similarly: 

r ^ t /" /(Ode 
J2 < const 



/ 



f| . vr ^||€||>i/||x-i/|| 

dr 



The final integral is finite if and only if a + f3 > d; and in this case, we have 
the estimate < J2 < const ■ \\x — y\\ a+ ^~ d . 
Finally, 

(4.33) ^<ll^ll<Vll^j 

2 /-Vlk-wll dr 
~\\ x -y\\ - J r a+P-d-l- 

We evaluate the integrals (4.32) and (4.33) for the different cases of a + j3 
to obtain the result. □ 

The following is an immediate consequence of Proposition 4.4 and The- 
orem 4.5. 

Corollary 4.7. Every random-field solution u given by Theorem 4-5 
has a continuous modification for all t > 0. 

We devote the remainder of this section to a special case of (4.1) namely 

(4.34) -u t (x) = (Au t )(x) + F t (x), 

where uq = 0, x G R 3 , t > 0, and the Laplacian acts on the x variable only. 
The noise F is a centered Gaussian noise, as before, that is white in time 
and homogeneous in space with a correlation function / that satisfies the 
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(4.37) "" Xll>1/M 



following for a fixed q 6 R: 

(4 ' 35> /K)S M» fo ^"KII>e. 

According to Theorem 3.17 on page 42, such correlation functions exist. 

The following lemma will be useful for the proof of the main result of 
this section. 

Lemma 4.8. If g : R 3 i-> R + is a Borel-measurable radial function, then 

(4.36) / (1 — cos(x • y))g{x) dx > const • / g(x)dx, 

■/|MI>i/ll»ll ^IWI>VIMI 

uniformly for all y G R 3 \ {0}. 
Proof. Clearly, 

/ (1 — cos(x • y))g{x) dx 

J\\x\\>l/\\y\\ 

/*oo /* 

= / r 2 R(r) dr / d9 (1 - cos(y • r0)), 

Jl/\\y\\ Js 2 

where R is the function on R + defined by i?(||x||) := g(x) for all x G R 3 . 
But for all r > 0, the d#-integral can be computed as 

(4.38) f (1 - cos(y • r9)) d6 = const • (l - Sm( ,^ ll) ^ , 

and this is bounded below uniformly, as long as r > l/||y||. We combine the 
preceding two displays to obtain the result. □ 

The following is the main result concerning (4.34). 

Theorem 4.9. Consider the stochastic heat equation (4.34) mR 3 , where 
the correlation function f of the noise satisfies (4.35) for a given fixed value 
ggR. Then: 

(a) (4.34) has a random-field solution u iff q > 1; 

(b) x i->- u t (x) has a continuous modification for all t > iff q > 2. 

Remark 4.10. Theorem 4.9, and general facts about stationary Gauss- 
ian processes [see Belyaev's dichotomy [MR06, Theorem 5.3.10, p. 213], for 
instance], together prove that when q £ (1,2], the stochastic heat equation 
(4.35) has a random-field solution u that almost surely has infinite oscilla- 
tions in every open space-time set. This example was mentioned at the end 
of Introduction. □ 
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Proof. In order to show existence of a random-field solution, it suffices 
to show that T(l) < oo if and only if q > 1. Because = ||£|| 2 , we may 
write T(l) as follows: 



T(l) = -L f - 

1 ; 8tt 3 y R3 1 



no 



(4.39) 87r3 7 R3 l + 2||e|| 2 

/i +h 

8vr 3 ' 

where 

/(o .,. ^ r ._ /• ho 



ncii>ei + 2|ieii 2 



(4-40) /i:= / /ffi dg and / 2 := / 

Direct inspection reveals that 

(4.41) h x / - r . dr and J 2 x / 

7o (logr)« A r(logr)9 

It follows readily from this that T(l) < 1 if and only if q > 1. 

We now turn our attention to the second part of the proof. Throughout, 
we assume that \\x — y\\ < 1/e, and consider the following integral: 



dr. 



(4.42) 



where 



8vr 3 



t ._ f 1 - cos(g • (x - y)) f 

Jl -4n<e i+2iieip /(0de ' 

(4.43) *:=/ lz^M/ (0 ^ 

J3:= / 1 Tnm f(0^. 

■/e<||£||<l/||x-j/|| i + 

We estimate each of the integral separately. The first term can be dealt with 
easily, and we obtain the following, using similar computations to those in 
the proof of Theorem 4.5: 

(4.44) Ji x \\x -y\\ 2 . 
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The estimation of the second term requires a little bit more work, viz., 
J,<U l-cos|.(x- y) ) /K)d; 
2 JU\\>l/\\x-y\\ llsll 



(4.45) 



- const • / d ^ 



lien 2 

< const • / - — — dr 



y\\ 

OO J 



'i/llx-i/ll r(logr)9 
const • I log - 



; \ -9+1 



|x - y\ 

Finally, we consider the final term J3: 



h~\\x-y\\ 2 [ /(0d£ 
(4.46) x || x _ y ||2^ 
< const • I log • 



" Ik-all r 

dr 

1 ^ 9 



\x - y\ 

Upon combining the above estimates we obtain the bound 



1 



-9+1 



(4.47) \d(x,y)\ 2 < const • log 

V \\ x ~y\\ 

Next, we compute a similar lower bound for \d(x,y)\ 2 . Since the inte- 
grands are nonnegative throughout, we may consider only J 2 . In that case, 
Lemma 4.8 yields the following: 

J 2 > const • / d£ 

(4.48) U 

> const • I log • 



no 

>l/\\x-y\\ 



\x - y\ 

Thus far, we have proved that 

(4.49) d(x,y)~\log(\\x-y\\)\ (1 -^ 2 , 
uniformly, as long as \\x — y\\ < 1/e. From this, we obtain 

(4.50) logN d (e)~e 2 ^-i\ 

valid for < e < 1/e. In particular, the metric-entropy condition (4.21) 
applies if and only if q > 2. Since the other conditions of Proposition 4.4 
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hold [for elementary reasons], the second part of the theorem follows from 
Proposition 4.4. □ 



CHAPTER 5 



The Nonlinear Equation 



The primary goal of this chapter is to study the fully- nonlinear stochastic 
heat equation (1.1) as described in the introduction. 

In the first part, we derive a series of a priori estimates that ultimately 
lead to the proof of Theorem 1.2. The latter theorem shows that the finite- 
potential Condition 1.1 is sufficient for the existence of a mild solution to 
the stochastic heat equation. As a byproduct, that theorem also yields a 
temporal growth rate for the solution. This means that under some natural 
conditions on the multiplicative nonlinearity a, the mild solution will not 
be intermittent. 

The second part is devoted to the proofs of Theorems 1.8 and 1.11, and 
thereby establishing the fact that, in contrast to the preceding discussion, if 
"there is enough symmetry and nonlinearity," then the mild solution to the 
stochastic heat equation is weakly intermittent. 

In the third and final part, we give a partial answer to a deep question 
of David Nualart who asked about the "effect of drift" on the intermittence 
of the solution. In particular, we show that if the drift is exactly linear — 
which corresponds to a massive and/or dissipative version of (1.1) — then 
there is frequently an explicit phase transition which describes the amount 
of drift needed in order to offset the intermittent multiplicative effect of the 
underlying noise. 



The main goal of this section is to prove Theorem 1.2. With that aim 
in mind, we can formulate (1.1), in mild form, as follows: 



5.1. Existence and Uniqueness 



(5.1) 
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As is customary, we seek to find a mild solution that satisfies the following 
integrability condition: 

(5.2) sup sup E (\ut(x)\ 2 ) < oo for all T > 0. 

te[0,T] x€R d 

In order to prove Theorem 1.2 we apply a familiar fixed-point argument, 
though the details of this argument are not entirely standard. 

Let J- := {J~t}t>o denote the right-continuous complete filtration gener- 
ated by the noise F. Specifically, for every positive t, we define J-J 1 to be the 
(T-algebra generated by random variables of the form of the Wiener integral 
Jj ot j xRd 4>s(x) F(dsdx), as <j) ranges over L 2 ([0,i] x H d ). Define T\ to be 
the P-completion of J 7 ®, and finally define 

(5.3) Tt'=[\T\ 

s>t 

as the right-continuous extension. 

We recall from Walsh [Wal86] that a random field {vt(x)} t>0 x€ Rd is 
predictable if it can be realized as an L 2 (P)-limit of finite linear combination 
of random fields of the type 

(5.4) zt(x)(uj) := X(u})l( a $ xA (t,x) for t > 0, x G R d , and u G f2, 

where: < a < b < oo, A C R d is compact, and X is an J^-measurable and 
bounded random variable. 1 Define for all predictable random fields v, 

(5.5) (Av) t (x):= p t - s {y-x)<j{v s (y))F(dsdy), 

J[0,t]x~R d 

and 

(5.6) (Bv) t (x) := ds dy p t - s (y - x)b(v s (y)), 

Jo Jn d 

provided that the integrals exist: The first integral must exist in the sense 
of Walsh [Wal86]; and the second in the sense of Lebesgue. 
Define for all /3,p > 0, and all predictable random fields v, 

(5.7) M\f3, P ■= sup sup 

t>o xeR d 

It is easy to see that the preceding defines a [pseudo-] norm on random 
fields, for every fixed choice of > 0. In fact, these are one among 



e-^E(\v t (x)n 



We are using the standard "(f2 , T , P)" notation of probability for the underlying proba- 
bility space, of course. 
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many possible infinite-dimensional L p -norms. And the corresponding L p - 
type space is denoted by Bp :P . We make the following definition which will 
be in force throughout the rest of the paper. 

Definition 5.1. Let Bp :P denote the collection of all [equivalence classes 
of modifications of] predictable random fields X := {X t (x)} t>0 xen d such 
that ||X||0 jP < oo. 

One can easily checks easily that || • \\p jP defines a pseudo-norm on Bp tP . 
Moreover, if we identify X G Bp jP with Y G Bp jP when \\X — Y\\p :P = 0, then 
[the resulting collection of equivalence classes in] Bp :P becomes a Banach 
space. Because \\X — Y\\p p = if and only if X and Y are modifications 
of one another, it follows that — after the usual identification of a process 
with its modifications — Bp iP is a Banach space of [equivalence classes of] 
functions with finite || • \\p p norm. 

Our next two lemmas contain a priori estimates on Walsh-type stochas- 
tic integrals, as well as certain Lebesgue integrals. Among other things, 
these lemmas show that B and A are bounded linear maps from predictable 
processes to predictable processes. These lemmas are motivated strongly 
by the theory of optimal regularity for parabolic equations, as is our entire 
approach to the proof of Theorem 1.2; see Lunardi [Lun95]. We follow the 
main idea of optimal regularity, and aim to find a good function space such 
that if uq resides in that function space, then ut has to live in the same 
function space for all t. As we shall soon see, the previously-defined Ba- 
nach spaces {Bp jP }p jP>0 form excellent candidates for those function spaces. 
In a rather different context, this general idea appears also in Dalang and 
Mueller [DM03]. Those authors show that L 2 (R d ) is also a good candidate 
for such a function space provided that cr(0) = 0. 

Here and throughout, we will use the following notation on Lipschitz 
functions. 

Convention 5.2. If g : TL d — > R is Lipschitz continuous, then we can 
find finite constants C g and D 9 such that 

(5.8) \g(x)\ <C g + D g \x\ for all x G R°\ 

To be concrete, we choose C g := \g(0)\ and D g := Lip fl , to be concrete. 

As mentioned above, the next two results describe a priori estimates for 
the Walsh-integral-processes Bv and Av when v is a nice predictable random 
field. Together, they imply that the random linear operators A and B map 
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each and every Bp lP into itself boundedly and continuously. The respective 
operator norms are both described in terms of a replica potential of the 
correlation function /. 



Lemma 5.3. For all integers p > 2, real numbers (3 > 0, and predictable 
random fields v and w, 

(5-9) IW^<| (^+DfclMk P ), 

and 

(5.10) \\Bv - Bw\\p tP < ^j^-\\v - w\\/3 tP . 

Proof. On one hand, the triangle inequality implies that E(\(Bu) t (x)\ p ) 
is bounded above by the following quantity: 

I dsi f dyx--- f ds p I dy p JJ p t - Sk (tk - x) ■ E [ JJ \b(v s {yj))\ J . 
Jo Jr* Jo Jn* ^ \f = \ J 

On the other hand, the generalized Holder inequality tells us that 

(5.H) E[ni6K.( % ))i) <nii6K.(%))n P . 

V=l / 3=1 

Therefore, we can conclude that 

(5.12) E(\(Bv) t (x)\ p )< (c b t + n b J*dsJ d yPt - 8 (y-x)\\v 8 (y)\\ p 

We multiply the preceding by exp(— fit) and take the the (l/p)-th root to 
find that 

i il v 

e-^E(\(Bv) t (x)\ p )\ 

< C b te-Wv + D 6 / ds dy Q-^-'^pt-siy ~ x)e-^\\v s {y)\\ P 
Jo Jn d 



(5.13) 



<c b te-^ + ^\\ v y, p . 

The first display of the lemma follows because texp(—(3t/p) < p/(ef3) for all 
t > 0. In order to obtain the second display we note that 

(5.14) \(Bv) t (x) - (Bw) t (x)\ < Lip, • (B^v - w\)) t (x), 
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where B\ is defined exactly as B was, but with b(x) replaced by b\ (x) = x. 
Because we may choose = and = 1, the second assertion of the 
lemma follows from the first. □ 

Lemma 5.4. For all even integers p > 2, real numbers (3 > 0, and 
predictable random fields v and w, 



(5.15) \\Av\\ PtP < z p (C CT + B a \\v\y tP ) y/(Bvp/ P f)(0), 
and 

(5.16) \\Av - Aw\\p tP < z p Lvp a \\v - w\\p <p yj (-R 2 /3/ P /)(0)- 

Proof. According to Davis's formulation [Dav76] of the Burkholder- 
Davis-Gundy inequality, E(\( K Av)t{x)\ p ) is bounded above by Zp times the 
expectation of 

ft r r p l 2 

(5.17) / ds / dy I dz V s (y , z)p t - s (y - x)p t - s (z - x)f(z - y) 
Jo Jn d JR d 

where 

(5.18) V s (y,z) := a(v s (y))a(v s (z)). 
By the generalized Holder inequality, 

/ P/2 \ p/2 

(5.19) E J[V aj {3,j,Zj) < [J IHvs^y^Wjaivs^z^Wp. 

\J=1 / 3=1 

Consequently, E (| (Av)t(x)\ p ) is bounded above by Zp times 

(5.20) ( [ ds [ dy f dzVt(y,z)pt-,(y-x)pt-,(z-x)f(z-y)) ' , 
\Jo JR d JR d J 

where 

(5-21) V^y ,z) := \\a(v s (y))\\p\WMz))\\ p . 

We can note that for all s > and y € H d , 

\\<T(v s (y))\\ p <C a + D a \\v s (y)\\ p 

(5.22) <C a + B a e^ p \\v\\^p 

<<J>'/p ( Cct + D <T ||«|| j8>p ). 

Therefore, a line or two of computation yield 

/ roo \p/2 

(5.23) E(\(Av) t (x)\ p )<e^z p (C a + B a \\v\y : p) p n e- 2 ^ p n s dsj , 
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where 



/R d Jn d 

(5 - 24) = (p;p./) (o) 



Hs := / da db f(a - b)p s (a)p s (b) 
Jn d Jn d 

{P*,Psf) (' 
(P/)(0). 



And hence, 

(5.25) E (|(^M*)n < (C CT + T> c \\v\\^f ((i? 2/3/p /)(0)) p/2 , 

The first assertion of the lemma follows immediately from this. 
In order to deduce the second assertion we note that 

(5.26) \\Av - Aw\\fo < Lip a • Ai(\v - w\), 

where Ai is the same as A, but with a(x) replaced by g\(x) = x. Therefore, 
the first assertion of the lemma implies the second. This completes the 
proof. □ 

We are ready to begin a more-or-less standard iterative construction that 
is used to prove Theorem 1.2. 
Let 

(5.27) u° t (x) := u (x), 
and define iteratively: For all n > 0, 

(5.28) u^\x) = (P t u )(x) + (Bu n ) t (x) + (Au n ) t (x). 
Lemma 5.5. Choose and fix ft > and an even integer p > 2. // 

(5.29) pRa + ZpB(T ^ { R 2f3/pfm <i, 

then sup n>0 ^"H^p < oo. 

Proof. Because Ptuo is bounded, uniformly in modulus, by sup xeR d \uq(x 
the triangle inequality implies that 

(5.30) \\u n+1 \\(3, p < sup |«o(x)| + \\Bu n \\ fi>p + \\Au n \\p, p . 

x£R d 

Lemmas 5.3 and 5.4, and a few lines of direct computation, together imply 
that 



(5.31) IK +1 |k P < A + B 
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where 

(5.32) A := sup \u (x)\ + ^ + z p C a J(R^ /p f)(0), 
and 

(5.33) B:=^ + z p B a ^(R 2f3/p f)(0). 
Iteration yields the bound 



(5.34) ||u n+1 Wfa < A 1 + B + ■ ■ ■ + B n ~ x + B n sup |u (. 



x 



X 



Consequently, if B < 1 then 

(5.35) supll^H^p < -. 

k>i i — a 

Since ||it ||/3 iP < sup xGR d |tto(^)| < °o> the lemma follows. □ 
We now have all the technical estimates for the proof of Theorem 1.2. 

5.1.1. Proof of Theorem 1.2. Without loss of generality, we can find 
(3 > such that Q(p,f3) < 1, otherwise there is nothing to prove. Choose 
and fix such a j3. 

Thanks to Lemma 5.5, every u n is well defined and ||ii n ||/3, p is finite, 
uniformly in n. In particular, u n G S/3, p - Next we apply Lemmas 5.3 and 
5.4 — with C g := |g(0)| and D g := Lip g — to find that 

\\u n+l - u n \\^ p < \\Bu n - Bu n -% iP + \\Au n - Au n -%, p 

(5.36) 

< ||u n -u n -%, P -Q(p,/3). 

Because Q(p,f3) < 1, the preceding implies that 



(5.37) - u"!!^ < const • J^{Q(p,/3)} n < oo. 

n=l n=l 

Therefore, we can find a predictable random field u°° G Sg, p such that 

(5.38) lim o ||u n -u oo || /3i p = 0. 
Our arguments can be adjusted to show also that 

(5.39) lim \\Bu n -Bu oo \\ 0p = O, 
as well as 

(5.40) lim \\Au n -Au°° | |/j, p = 0. 
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This proves that u°° is another solution to (5.1). As we have mentioned, u 
is the almost-surely unique solution to (5.1). Therefore, u°° is equal to u, up 
to evanescence. It follows that u G Bp, P for all (3 > such that Q(p,f3) < 1. 
This proves the theorem. □ 



5.2. Intermittency 

In this section we prove Theorems 1.8 and 1.11, which state that the 
solution to the stochastic heat equation can be weakly intermittent in the 
presence of enough symmetry and nonlinearity. The basic idea is to follow 
our earlier work on space-time white noise [FK09] and apply the renewal- 
theory methods of Choquet and Deny [CD60]. However, this turns out to 
be a difficult adaptation which appears to require a fair bit of harmonic 
analysis. 

It might help to recall the definition (1.15) [page 4] of the function T, 
and the relation [Theorem 1.5, page 5] between the positive number T(/3) 
and the maximum value of the replica /3-potential of the correlation function 
/, as well as the positive number (^ ( g/)(0) of that replica potential at the 
origin. 

Proof of Theorem 1.8. We can assume, without loss of generality, 
that there exists /3 > such that 

(5-41) T(/3) > A 

for there is nothing left to prove otherwise. Now consider any such (3. 

Since b(x) = 0, no > rj, and a(u) > L a \u\, we can apply (5.1) to deduce 
that for all x, y G H d and t > 0, 

E(\u t (x)u t (y)\) 

(5.42) 

>E(u t (x)u t (y)) 

> v 2 +[ds[ dz f dz' W s (z,z')p t - s (z-x)pts(z' -y)f(z-z'), 
JO JR d JR d 

where 

W s (z,z') :=E(a(u s (z))a(u s (z')) 
(5 ' 43) >L 2 a E (\u s (z)u s (z')\). 
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Consider the following R + - valued functions on (R d ) 2 : 

poo 

H p [a,b) := / e-^E(\u t (a)u t (b)\)dt, 
Jo 

(5 - 44) Gp(a,b):= re-^p t (a)p t (b)dt, 

Jo 

F(a,b) :=f(a-b). 

Also consider the linear operator Ap defined as follows: For all nonnegative 
Bor el-measurable functions h : (R d ) 2 — > R + , 

(5.45) (Aph)(x , y) := (fIi * Gp) (x , y). 

A line or two of computation shows that the preceding is simply a quick way 
to write the following: 

(Aph){x,y) 

(5.46) = j Rd da J Rd db F ( a MK a M G p( x - a ,y - b) 

= [ e~^dt [ da [ db f(a - b)h(a ,b)p t (x - a)p t (y - a). 
Jo JR d Jn d 

With the preceding definitions under way, we can write equation (5.42) 

in short hand as follows: 

(5.47) H p {x,y)>^ + Ll (ApHp) (x,y). 

Since F > on (R d ) 2 , we can apply the preceding to find the following 
pointwise bounds: 

Hp > £ 1 + Ll (Ap + Ll(ApH^ 
= V jl + L^jApl + LiA}Hp, 
> V jl + Ll V jApl + LiAp- (^+LlApHp 
= |l + L^i,l + L^^l + L^, 



(5.48) 
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where l(x,y) := 1 for all x,y € H d . By applying induction we may arrive 
at the following simplified bound: 

(5.49) £=0 

Because the parameter k > is arbitrary, it follows that 

2 00 

(5.50) Hp^-.^U^l. 

p e=o 

In order to better understand the behavior of this infinite sum, we begin by 
inspecting only the first few terms. 

The first term in the sum is identically 1. And the more interesting 
second term can be written as \? a multipled by 

(Apl)(x,y) = (F*G fi ) (x,y) 

f'OC P p 

(5 51) =/ e~ m dt da db f(a - b)p t (a - x)p t (b - y) 

x ' JO JR d JR d 

> j±p J" e~*dt J dt /(Oe-^Oe*^), 

thanks to Proposition 3.6 [p. 31]. We change the order of the double integral 
[dt d£] to find that 

(5.52) {A , 1)(x , v) >—j—m. i( . 

[Theorem 1.5 and condition (1.1) together justify the use of Fubini's theo- 
rem.] 

In order to bound the third term in the infinite sum in (5.50), we need 
to estimate the following quantity: 

(A}l)(x,y) = (F(Apl)*Gp) (x,y) 

(5 - 53) ^ i f jmz 7( , 

- (2vr) d J (3 + 2Re*(0 ^ 

where 

pOO p P 

(5.54) Z^(x,y) := J e'^dt J da J db f \a - b)S< a S t {a ~ x)pt(b - y) . 
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[The same ideas that were applied to the second term can be applied here, 
in exactly the same manner, to produce this bound.] The Fourier transform 
of the function 

(5.55) a i— > exp(i£ • a)pt(a — x) 
is 

(5.56) C^exp(^ + C)-x-^ + 0). 
Therefore, Proposition 3.6 implies the following: 

Z 6 (x,y) e-% J dC /( C )e«(€-K)-(x-y)-2ffla*(€+0 

= J (3 + 2Re*(e + C) dC ' 



And therefore 

12 



,»(£i+62)-(x-j/) 



(2tt) m 7 /? + 2Re^(£i) J /3 + 2Re*(^+6) 
Now, we can apply induction to deduce that we have the following esti- 
mate [used to analyse the £th term in the infinite sum in (5.50)] in general: 
For all integers t > 1 and x, y £ R rf , 

> _j_ / /(a)da / /(6)d^ 2 

(5.59) - ( 27r )« J/3 + 2Re^(6) 7 /3 + 2Retf (& + &) ' 

A&W p <Ej = iO-(»-y) 



/3 + 2Rett(& + ■■■+&) 
Although the preceding multiple integral is manifestly nonnegative, the 
integrand itself is complex valued. Fortunately, we wish to only understand 
the behavior of F^l on the diagonal of (R d ) 2 . In that case, the integrand is 
real and nonnegative. As such, we can estimate the integrand directly. In 
order to do so, let us set y := x in (5.59), and then plug the result in (5.50), 
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to find that 

poo 

inf / e-^E(|ut(a;)| 2 )dt 

x£~R d J 

(5.60) 

^^•E(4^y Rd d ^---y Rd d ^ n (/3+2Re ^/ + ... +0)) ' 

A change of variables yields the following: 



poo 

inf / e-^E(\u t (x)\ 2 )dt 
eR d Jo 



xeR d Jo 

( 5 - 61 ) v 2 oo ^ < ^, 

" ^ • g L ^ " ' L ^ 09T 2Re*(^)) ' 

where zo := 0. 

The preceding holds even without Condition 1.7 [p. 6]. But now we 
recall that Condition 1.7 is in place, and use it to produce the announced 
lower bound on inf x£K d 7 X (2). 

Define 

(5.62) S := |x := (xi , . . . , x d ) G R rf : sign(xi) = • • • = sign(x rf )| . 

Equivalently, £ := R/[ U Rl. Because the terms under the product sign 
in (5.61) are all individually nonnegative, Condition 1.7 assures us that the 
following holds: 



poo 

inf / e-^E(|«t(x)| 2 )dt 
eR d Jo 



xeR d Jo 

If 21, . . . , ^ € S, then the absolute value of the fcth coordinate of Zj — Zj-i 
is less than or equal to the absolute value of the kth coordinate of Zj for all 
k = 1 , . . . d; therefore 

(5.64) f( Zj - zj-t) > f{zj), 
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thanks to Condition 1.7. Consequently, 

POO 

inf / e-^E(\u t (x)\ 2 )dt 
x&n d Jo 



(5-65) " /? ^(27r) M 7 E , a l = l(^ + 2Re*( 



In particular, if there exists /3 > such that 
(5 66) _L /' Mjg > ] -2 

then 

(5.67) J e _/3 *E (\u t (x)\ 2 ^ dt = oo for all x G H d . 

Thanks to symmetry considerations, Condition 1.7 has the following conse- 
quence: 

f(0^ n - d+ i f /(0d£ 



7 S /3 + 2Re*(£) u 



(5.68) i s /3 + 2Re*(0 ~ V /3 + 2Re* (0 

= 2- d+1 T{p). 

This and (5.66) together imply that (5.67) holds whenever 

2<z-i 

(5.69) T(/3) > 



Now we can apply a real-variable argument to prove that if (5.67) holds 
for some (3 > 0, then 

(5.70) 7 X (2) > /3 for all x G R d . 

This and (5.41) together imply the theorem. 

Indeed, let us suppose to the contrary that (5.67) holds for our (3, and 
yet mf x&R d 7 X (2) < j3 for the very same j3. It follows immediately that there 
exists x G R d , <5 G (0,13), and C G (0 , oo) such that 

(5.71) E (\u t (x)\ 2 ) < Ce ((3 ~ 5)t for all t > 0. 

And hence, (5.67) cannot hold in this case. This produces a contradiction, 
and shows that (5.67) implies the theorem. □ 
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Proof of Theorem 1.11. We begin as we did with (5.42), but can no 
longer apply the inequality in (5.43). To circumvent that, note that for all 
qo € (0,(?) there exists A := A(qo) € (0 , oo) such that a(y) > qo\y\ as soon 
as \y\ > A. We have assumed that P{u s (y) > 0} = 1 for all s > and 
y £~R d . Therefore, 

W s (z , z') > q%E (u s (z)u s (z') ; u s {z) A u s {z) > A) 

> q 2 E (u s (z)u s (z')) - q 2 A 2 - q AE(u s (z) ; u s (z) > A) 

(5.72) -q AE(u s (z') ; u s (z') > A) 

> qlE (u s (z)u s (z')) - qlA 2 - q A {E (u s (z)) + E (u s (z')) } . 

On the other hand, (5.1) guarantees that 

0<E(u t (x)) 

(5.73) = (P t «o) (x) 

< ||«o||l°°(r<*)- 

Consequently, 

(5.74) W s (z , z') > ql {E (u s (z)u s (z')) -A,}, 
where 

(5.75) A* := max ^A 2 , 2||uo|| I/ oo ( - R d^ . 

Now we apply the recursion argument in the proof of Theorem 1.8, and find 
the following pointwise bounds: 

Hp > j + ql {ApHp - A.Apl} 

2 / 2 \ 

>j + Q 2 o(j~ Apl + q$A%Hfi - q 4 Q A*A}l 

(5.76) 



2 / 2 \ N 

valid for every integer N > 1. We apply the following obvious inequalities: 
rj 2 / (3 > rj 2 / (3 — A* to the first term on the right; and Hp > rj 2 / '(3 to the last 
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bracketed term, to find that for all integers N > 1, 

/ 2 \ N+ 1 

(5-77) H,>{ T L-AAY j ^^1. 

We can now let N f oo and apply the same estimate that we used to derive 
(5.65), and deduce the following bound: 



poo 

inf / e- /3t E(|u i (x)| 2 )dt 
xen d Jo 



(5.78) 



I V 2 \ °° I Qo f f( z ) 
~ ' ~ V " S {WF L (0 + 2Re^(z)) d " y 



' 2 ^E^^))' 



£=0 

Thanks to Theorem 1.5, the preceding implies the following bound: 

poo / 2 \ 00 / 2 \ J 

(5.79) mf^ e"*E (M*)| 2 ) dt > - • £ (^T W)(0)J 

Because (Rof )(0) = oo, we can find a /?o > such that 

2 d-i 



(5.80) (%,/)(0) > 2 • 
For that choice of /3o, the preceding sum diverges. Therefore, 

;>oo 

(5.81) / e~^ ot E(\u t (x)\ 2 ) dt = oo provided that r] > ^/J^Xl. 
Jo 

This and an elementary real-variable argument; (5.71) and together prove 
the theorem. □ 



5.3. The Massive and Dissipative Operators 

David Nualart asked us about the effect of the drift coefficient b in (1.1) 
on the intermittent behavior of the solution to the stochastic heat equation 
(1.1). At present, we have only an answer to this in a special but physically- 
interesting family of cases. 

Indeed, let us consider the stochastic heat equation 

(5.82) ^u t (x) = (Cu t )(x) + K t {x) + a(u t (x))F t (x), 
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where x G R d , t > 0, A G R, and F is as before. Moreover, a : R — > R is 
Lipschitz continuous, also as before. That is, (5.82) corresponds to the drift- 
free stochastic heat equation for the massive operator CW := C + (A/2)/ 
when A > 0, the dissipative operator = £ — \\/2\I when A < 0, and 
the free operator = C when A = 0. 2 Of course, Theorem 1.2 — applied 
with b(u) := \u/2 — guarantees us of the existence and uniqueness of a mild 
solution to (5.82). 

The operator is the generator of the semigroup {-P/^}t>o defined 

by 

(5.83) {Pi X) <t>){x):=e xt '\P t ct>){x). 

This can be seen immediately by a semi-formal differential of t \— > at t = 
0; and it is easy to make the argument rigorous as well. The corresponding 
"transition functions" are given by 

(5.84) P? ) (v-x):=e xt ' 2 p t (y-x). 

Then the domain Dom[£( A '] of the definition of is the same as 
Dom[£], and 

(5.85) £ (A) = £<p + ^4> for all <j> G Dom[£ (A) ]. 

Let {Pt^}t>o denote the adjoint [or dual, in probabilistic terms] semi- 
group. That is, 

(5.86) (P; {A V)W:=e At / 2 (P^)(x). 
with corresponding transition functions, 

(5.87) p* t iX) (y - x) := e u ' 2 p t {x - y). 

And finally there is also a corresponding replica semigroup, 

(5.88) (lfV)(*):=e At (P^)(*), 
whose resolvent is described by the following: 

POO 

(5.89) {R {x) cj)){x) := e-( a - x > (P s cj))(x) ds for all a > A. 

Jo 

We might note that f = R a -\f is merely a shift of the the free replica 
resolvent of /. Therefore, the proof of Theorem 1.2 goes through unhindered, 
and after accounting for the mentioned shift, produces the following: 

2 The operator jC^ is also known as the "relativistic" form of C. 
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Theorem 5.6. Suppose uq : H d — > R is bounded and measurable. Then, 
under Condition 1.1, the mild solution to (5.82) satisfies the following: For 
all integers p > 2 and A 6 R, 

(5.90) %(p) < A + V - inf J a > : (R a f)(0) < ' " 



where z p denotes the largest positive zero of the Hermite polynomial He p . 
Recall the function Q of Theorem 1.2. Since 

(5.91) Q(p,0) > max , z p Lip a ^ (R 2 ? /p f)(0)) 

a few lines of arithmetic show that Theorem 5.6 provides us with a better 
upper bound than Theorem 1.2 for the top Liapounov L p -exponent of the 
mild solution to (1.1). Next, we produce instances where the solution is 
intermittent. 

First of all, note that according to (5.90), 

(5.92) %(2) < A + inf{a>0: {R a f){0) < 

because z 2 = 2. We apply similar "shifting arguments" together with Theo- 
rem 1.8 to deduce that the following offers a converse, under some symmetry 
and regularity conditions. We note, in advance, that when d = 1, the pre- 
ceding estimate and the following essentially match up. 

Theorem 5.7. Suppose that both Conditions 1.1 and 1.7 hold, n := 
inf^gj^d uo(x) > 0, and there exists L a £ (0 , oo) such that a(z) > L a \z\ for 
all z € R. Then, 

(5.93) inf 7x(2) > A + supL>0: (i? Q /)(0) > ^ 
where sup0 := 0. 

We end this chapter with the example mentioned on page 10 of the 
Introduction. 

Example 5.8. Consider the case that C = -(-A)«/ 2 is a power of the 
Laplacian. In that case, C is the generator of an isotropic stable process 
of index q, and q £ (0,2] necessarily. Consider also the case that f(x) = 
|| x ||-d+6 ^ a j>j esz kernel, where b G (0,d). Then [see (3.25) on page 30], 

(5,4) /(f) = where C^^^MV. 
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And therefore, 

(R a f)(0) = T(q) [by Theorem 1.5] 

C d>b [ U\\~ b 
(5.95) " (2vr)^ ,/ Rd a + 2||£||" ? 

C d , b a- 1+ ( d - b ^ f dz 



In other words, 



(2jr) d • 2( d - b )/i J Rd 



(5.96) (B«/)(0) 
where 



(5.97) V := — b - and 2l d(?fe := 7 °f • / 



d - b , „, Crf.t /" dz 

/ Rd Ikll 6 + IMI 9+6 ' 

Since 6 £ (0,d), 2l<2q& — and hence (R a f)(0) — is finite if and only if 
q + b > d; this is the sufficient condition for the existence of a unique mild 
solution to the resulting stochastic PDE. And not surprisingly, this "q + b > 
d" condition is the necessary and sufficient condition for the existence of 
a solution to the linear equation [Theorem 4.5, p. 48]. Moreover, when 
q + b > d, we can apply Theorems 5.6 and 5.7 to find that for all x 6 R d , 

(or t 2 \ v ) 
<7,(2)<A + (2l^Lip^) 1/(1 ^. 

In particular, consider the massive/dissipative "parabolic Anderson model," 

(5.99) = " (( A )^ /2u *) ( x ) + 2 Ut ^ + ««t( a; )^t( a; )> 

where no : R — > R is measurable and bounded uniformly away from zero 
and infinity, and ft / 0. The preceding discussion shows that the parabolic 
Anderson model has a solution if q + b > d. And when q + b > d, we obtain 
the following bounds for the upper L 2 -Liapounov exponent of the solution: 
For all x € R d , 

(or 2 \ ") 

^f^ 1 ) <7,(2)<A+(2l d ^ 2 ) 1/(1 ^. 

Theorem 1.2 shows also that 7*(p) < oo for all p > 2. Therefore, we 
have no weak intermittency if A < — (2l c z,<?,ftft 2 ) 1 / ( ' 1 ~ l/> , whereas there is weak 
intermittency if A > — (Std,^/* 2 /^ -1 )) 1 ^ 1-1 ^ . Our condition is sharp when, 
and only when, d = 1. In that one-dimensional case, we have a solution if 
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and only if q + b > 1, and if this inequality holds then 

= C lt y°° dz 

(5,01) ^-^7o * + ^ 

= ^ . B ( 1=* , 1 - 1=^ ), 
2"irq V i i J 

where B(x , y) := T(x)T(y)/T(x + y) is the beta function. Thanks to (5.94), 
this implies that 

< 5 - io2 > "^- ^-lU 8 

and we find that weak intermittency holds if and only if 
(5.103) A>-(^, 6 k 2 ) 1/M - 

Another simple though tedious computation shows that this example [ap- 
plied with q := 2] includes the material that led to (1.26) of page 10. 



CHAPTER 6 



Temperate Solutions to Parabolic SPDEs 

In this chapter we continue the study of the stochastic heat equation of 
the following form: 

(6.1) J^ t (x) = (Cu t (x)) + a(u t (x))F t (x), 

where a : H d — > R is Lipschitz continuous, as before. But here no can be a 
finite Borel measure. 

Such a problem arises naturally in various specific cases; see, for ex- 
ample, Bertini and Cancrini [BC95], Carmona and Molchanov [CM94], 
and Molchanov [Mol91]. It is also motivated by the physical literature on 
statistical mechanics, where uq typically represents the initial [probability] 
distribution of a particle system in a random environment that interacts 
with it random environment. In fact, the role of uq as a measure is typically 
taken for granted and little or no mention of uq is made explicitly in the 
physics literature [KS91, Kar87, KPZ86]. 

Thus, our present goal is to continue and produce solutions to (6.1), 
even though the initial data is a measure. 

The fact that uo is no longer a function poses some fundamental prob- 
lems for the existing theory of SPDEs. In fact, we argue next that in order 
to make sense of (6.1) in the case that uq is a finite Borel measure, we need 
to impose regularity conditions on both the semigroup {Pt}t>o and the ini- 
tial measure. More importantly, we need to develop a slightly less restrictive 
concept of a solution than the well known, as well as standard, notion of a 
mild solution. 

Let us begin with a few observations. 

It follows from the properties of the Walsh stochastic integral [Wal86, 
Chapter 2] that if it is a predictable random field which is the mild solution 
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of (5.1) [with 6 = 0] and satisfies (5.2), then 

e(m*)| 2 ) 

f ft f 

,(y - x)a(u s (y)) F(ds dy) 




> \(P t u )(x)\ 2 , 
where as usual, 

(6.3) (P t u )(x) := p t (y - x)u (dy). 

Thus, in particular, (5.2) implies that the following is a necessary condition 
for the existence of a mild solution: 

(6.4) sup sup \(P t u )(x)\ < oo for all T > 0. 
*6(0,T) x&~R d 

It is clear that the latter condition rules out many interesting choices of uq. 
For instance, consider the case that C = A and one of the most natural 
choices uq := 5 Z for initial data. In that case, 

-||a- Z || 2 /(4i) 

(6.5) (PA)(a) = -jjjps-, 
and hence 

(6.6) sup sup \(P t S z ){x)\ = oo for all T > and z G R rf . 

te(0,T) x£R d 

Thus, we can never have a mild solution to (6.1) when the initial data is 
a point mass. But when a(u) = const • u and C = A, (6.1) is expected to 
have a solution defined by the Feynman-Kac formula in many cases where 
/ is "nice"; see Bertini and Cancrini [BC95] and Carmona and Molchanov 
[CM94], for example. It is therefore natural to try to use a different defini- 
tion of a solution to handle such problems. 

It turns out that one can turn the suggested strategy into a fruitful 
rigorous approach. In fact, in this chapter we consider a different notion of 
solution which satisfies, among other things, the following condition in place 
of the more restrictive condition (5.2): 

(6.7) sup E (\ut(x)\ 2 ^ < oo for almost every t > 0. 

Now suppose we want to know when (6.1), with uq = 5 Z , has a random- 
field solution which satisfies (6.7). Similar considerations as those of the 
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previous paragraph tell us that a necessary condition is that pt is a bounded 
function for almost all t > 0. Thanks to Hawkes's theorem [Proposition 2.3, 
p. 19], (6.7) implies that 

(6.8) exp(-Retf) G L\K d ) for all t > 0. 

Thus, we will will need to assume, at the very least, that (6.8) holds. 

Recall that a mild solution u of (6.1) is a predictable random field that 
satisfies (5.1) with 6 = 0. To be specific, one requires that for all nonrandom 
pairs (t,x) G (0,oo) x R d , the random equation (5.1) holds almost surely. 
The following is a slightly less stringent notion of a solution. 

Definition 6.1. Let u := {ut(x)}t>o,xeK d be a predictable random 
field. We say that u is a temperate solution to (6.1) if there exists a null set 
iVo C R+ such that for all t g" No, the following holds for every x G Yi d : 

(6.9) u t (x) = (P t uo)(x) + / pt- s (y - x)a(u s (y))F(dsdy) a.s. 

Jo Jn d 

As we shall see in the next section, the stochastic integral in (6.9) can 
be defined properly though, technically speaking, it is not a Walsh integral. 
Regardless, we have the following, which is the main result of this chapter. 

Theorem 6.2. Suppose a : H d — > R is Lipschitz continuous and uq is 
a finite Borel measure on H d such that 

,0.10) / 'f° ( l>' + /«» g < co. 

1 ' Jw l + 2Re*({) 5 

//, in addition, (6.8) holds, then there exists a temperate solution u := 

{ut(x)} t >o,x€R, d t0 

There is also a corresponding [limited] uniqueness result for these tem- 
perate solutions; see Proposition 6.14. 

Let us remark that (6.10) holds if and only if: 

(a) Condition (1.1) holds; and 

(b) 



(6.11) f lit ° m 

1 ' y Rd l + 2Re*(0 



d£ < oo. 



In addition, in the case that uq is a positive-definite function, the preceding 
display holds if and only if uq has bounded potentials; more precisely, that 
(R a uo)(Q) = sup xe - R d(R a uo)(x) < oo for all a > 0; see Theorem 1.5 for a 
more precise statement. 
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Corollary 6.3. Suppose d = 1, a : R — > R is Lipschitz continuous, 
and Condition 1.1 and the following are both valid: 

da 



/oo 
T- 
-oo 1 



< OO. 



+ 2Rett(0 

Then, (6.1) has a temperate solution for every finite initial measure uq. 

Remark 6.4. Condition (6.12) is equivalent to the property that the 
replica process X has local times; see Hawkes [Haw86]. In that case, it can 
be shown, as in our earlier work with E. Nualart [FKN09], that (6.1) has a 
temperate solution even when F is space-time white noise. We have stated 
Corollary 6.3 for d = 1 only because (2.14) [see p. 15] implies fairly readily 
that (6.12) can never hold when d > 2. □ 

Remark 6.5. Corollary 6.3 has the following remarkable consequence: 
When (6.12) holds, the stochastic heat equation (6.1) has a temperate so- 
lution for every finite measure uq; this includes every uq £ L 1 (R). This 
property fails to hold for a mild solution that satisfies the usual integrability 
condition (5.2). □ 

We conclude this section with our final main result concerning temperate 
solutions to (6.1). This result yields an upper bound for the growth of the 
solution. Since our solution is not in general a mild one, we will need to also 
adapt our previously- used notion of Liapounov exponents. Thus, we intro- 
duce the following, which seems to have a number of desirable mathematical 
properties. 

Definition 6.6. We define the integrated L P (P) -Liapounov exponent 
\*(p) of the temperate solution to (6.1) — when it exists — as 

2/p 

P > : 



(6.13) \*(p):=mf{ 



where inf := oo. 



/ e _/M <^ sup E(\u t (x)\ p ) \ at <oo 
Jo [xen d J 



> ■ 



Theorem 6.7. Suppose (6.8) holds, a : R — > R is Lipschitz continu- 
ous and nonrandom, and uq is a nonrandom finite Borel measure on R d 
that satisfies (6.10). Then the temperate solution u to (6.1) satisfies the 
following: 

(6.14) A*(p)<infi/3>0: {Rpf)(0) < ' 



2z|Lip2 
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According to Theorem 1.5, the integral condition (6.10) implies that 
(R a f)(0) < oo for all a > 0, whence lim | g_ > . oo (^ j g/)(0) = 0, thanks to 
the dominated convergence theorem. Consequently, Theorem 6.7 implies, 
among other things, that \*(p) < oo for all p > 2. 

It would be interesting to find nontrivial conditions on / and C that 
ensure that the preceding temperate solution satisfies A* (2) > 0, and in this 
way derive a version of weak intermittency for temperate solutions. Such 
an undertaking would also require nondegeneracy conditions on no. But it 
seems hard to find nontrivial conditions that can be placed on an initial 
measure no that guarantee the strict positivity of A* (2) without assuming 
that no is a bounded measurable function [and not just a measure] which is 
also bounded away from zero. 

We end this section by making a final observation about weak intermit- 
tency. 

Suppose (Rof)(0) < oo. Then, A*(p) = provided that 
(6.15) Lip CT < ' 



^2z p 2(i? /)(0)' 



We can apply (6.15) with p = 2 to deduce that if a is sufficiently smooth in 
the sense that 

(6.16) Li Pa < 1 

V / 2( J Ro/)(0) 

then A* (2) = 0, and "weak intermittency" [with respect to the new Lia- 
pounov exponents {^*(p)} P >2\ does not hold. 



6.1. Stochastic Convolutions 

Define for all predictable random fields Z := {Zt{%)}t>o,xeR. d i 

(6.17) (p*ZF) t (x):= f [ p t - s {y-x)Z s (y)F(dsdy), 

JO JR d 

provided that the preceding Walsh-type stochastic integral exists. Note, in 
particular, that whenever v := {vt(x)} t>0 xeR d is a predictable random field, 

(6.18) Av = p * (cr o V )F, 

provided that the Walsh stochastic integrals are well defined. 
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According to the theory of Walsh [Wal86], the stochastic-convolution 
process p * ZF is well defined provided that the quantity 

(6.19) f ds f dy f dzE{\Z s (y)Z s (z)\)p^ s (y-x)p t _ s (z-x)f(y-z) 
JO JR d JR d 

is finite for all choices of (t , x) € (0 , oo) x R d . In that case, 
(p*ZF) t (x) 2 



E 
(6.20) 



= ds dy dzE(Z s (y)Z s (z))pt- s (y-x)p t - s (z-x)f(y-z). 
Jo Jn d Jn d 

Define for all > and all predictable random fields v := {t)((i)} (>0jl , gRl i, 



1/2 



(6.21) Mp{v) := ( f°° sup E (|^(x)| 2 ) dt 

\Jo xen d / 

Each Np is a norm on equivalence classes of square-integrable predictable 
processes that are modifications of one another. 



Definition 6.8. Let Z := {Z t (x)} t>0xeR d be a random field. We say 
that Z is (p, F)- Walsh integrable if Z is predictable and the quantity in 
(6.19) is finite. We frequently abuse notation and write "Walsh integrable" 
in place of "(p, F)-Walsh integrable." 



In other words, Walsh-integrable random fields are random fields for 
which the stochastic-convolution process p * ZF is defined by the stochastic 
integration theory of Walsh [Wal86] . 

The following is a key "embedding" theorem for our extension of sto- 
chastic convolutions. 



LEMMA 6.9. If Z := {Z t (x)} t>0 is Walsh integrable, then 
(6.22) Afp(p*ZF)<Afp(Z)-yJ(Rpf)(0) for all /3 > 0. 



E 
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Proof. Owing to (6.20), the following is valid for all £ > and x G R d : 

2 



(p*ZF) t (x) 



< [ sup E (\Z s (a)\ 2 ) ds [ dy [ dz p t _ s (y - x)p t _ s (z - x)f(y - z) 
Jo aeR d v 7 Jn d Jn d 

(6.23) 

= f sup E(\Z s (a)\ 2 ) (P t - s /)(0)da. 

Since the right-most quantity is independent of x G R d , we can take the 
supremum over all x G R ', multiply the preceding by exp(— fit) and then 
integrate [dt] to deduce the lemma. □ 



We now use Lemma 6.9 to extend the Walsh definition of a stochastic 
convolution as follows. Suppose (R a f)(0) < oo for some a > 0; because of 
Theorem 1.5 we know that (Rpf)(0) < oo for all [3 > 0. Define L 2 ^ to be the 
collection of all [equivalence classes of modifications of] predictable random 
fields Z such that Np{Z) < oo. 1 Clearly, L 2 ^ is a Banach space for every 
[3 > 0. 

Lemma 6.10. If Z e /or some /3 > 0, then there exist Z 1 , Z 2 , ■ ■ ■ e L 2 ^ 
such that lim ra „ i , 00 Z n = Z in L 2 ^, and Z n is Walsh integrable for all n > 1. 

Proof. We may assume without loss of generality that Z t (x) > for 
all £ > and x G R rf ; for otherwise we can consider the predictable random 
fields {Z t (x)) + and {Z t (x))~ separately [they too are in L 2 ]. 

For all t > 0, x G R d , and n > 1, let Z™(x) denote the minimum of ^(x) 
and n. Then Z ra converges in L\ to Z as n — >■ oo, thanks to the dominated 



^As is customary in the study of L p spaces, we treat the elements of slightly carelessly 
as if they are random fields, rather than equivalence classes of random fields. This is done, 
as in the case of L p spaces, merely to simplify the otherwise-cumbersome notation. There 
should be no loss in precision, as ought to be clear from the context. 
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convergence theorem. And 

fds! dy [ dzE(\Z™(y)Zl l (z)\)p t „ s (y-x)p t „ s (z-x)f(y-z) 
JO JR d JR d 

<n 2 ds dy dz p s (y - x)p s (z - x)f(y - z) 
Jo Jn d Jn d 

(6.24) 

= n 2 [ (P s f)(0)ds 
Jo 

< nV^/XO), 

which is finite. □ 

Now let us choose and fix some Z G Z/g. According to Lemma 6.10 
we can find Z n — >■ Z [in L 2 ] such that every Z n is (p, F)-Walsh integrable. 
Lemma 6.9 tells us that limn^oo (p * Z n F) exists in _L|. Thus, we can define 
the stochastic-convolution process p * ZF for every Z G L 2 as follows: 

(6.25) ip*ZF) t (x) := lim (p * Z n F) t (x), 

n— >oo 

where the limit takes place in Z/|. Of course, the preceding limit does not 
depend on the choice of {Z n }^ =1 . But it also does not depend on (5 because 

(6.26) L }^ L l if a < (3. 



We abuse notation slightly and write for all t > and x G R , 

(6.27) f f p t - s {y-x)Z s {y)F{dsdy) : = (p*ZF) t (x) for Z G M L 
Jo Jn d 



2 

0- 

/3>0 



It is easy to see that many of the standard properties of the Walsh stochastic 
convolution transfer, by limiting arguments, to those of the present exten- 
sion. Chief among them is the fact that Z i— > (p * ZF) is a random linear 
map from L 2 to itself. In other words, if Z, W G Lp and a, b G R, then 

(6.28) p * {aZ + bW)F = a{p * ZF) + b(p * ZF). 



We emphasize that the preceding means that, as elements of Z/|, the two 
sides agree. In particular, because of the Fubini theorem, we obtain the 
following equivalent formulation of (6.28): There exists a null set No C R+ 
such that for all t N and all x G R d , 

(6.29) P { (p * (aZ + bW) Fj ^ (x) = a (p * ZF^j ^x) + b(p* WF^j ^ (x) } 
is equal to one. 
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We conclude this section by emphasizing that the stochastic convolution 
process i H> (p * ZF) t is hereby not defined for all t > 0. However, it is 
defined for almost all t > 0; and if (p * ZF) t is defined for a given t > 0, 
then {p * ZF)t(x) is well-defined for every x 6 R d , and x i->- (p * ZF) t (x) is 
a random field in the usual sense. 

Our stochastic convolution is thus quite different from the existing ones 
in the literature. For example, for the most commonly-used infinite-dimensional 
stochastic convolution [see, for example, Chapter 5 of Da Prato and Zabczyk 
[DPZ92, §5.1.2]], each random process x H> (p*ZF) t (x) is defined for every 
t > 0, typically as an element of a certain Hilbert space. But, in that case, 
the random variables (p* ZF) t (x) cannot in general be defined for every x. 



6.2. A Priori Estimates 

We now proceed to establish a priori estimates for Walsh-type stochastic 
integrals of the form Av, thereby showing that the random linear operator A 
maps every L 2 ^ to itself continuously and boundedly. As such, the following 
result should be compared to Lemma 5.4 on page 59. 

Lemma 6.11. For all (3 > and predictable random fields v and w, 
(6.30) M p {Av) < + BvA/>(^ v^TKO), 

and 



(6.31) A/> {Av - Aw) < Lip CT A/>( V - w) y / 2(^/)(0). 

Proof. According to (5.20) and (5.22), E(\(Av) t (x)\ 2 ) is bounded above 

by 

Q 2 s ds dy dz p t - s (y - x)p t _ s (z - x)f(z - y) 

(6.32) = folds f dy I dzp t - s (y)p t - s (z)f(z-y) 

Jo Jn d Jn d 

(6.33) = f Q 2 s -(P t - s f)(0)ds, 

Jo 

where 

(6.34) Q s :=C a + D a sup K(a)|| 2 . 

aeR d 
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Therefore 

poo 



(6.35) 



Since 



poo /*oo roo 

/ e-^E(|(^) t (x)| 2 ) dt< / e-P s Q 2 s &s- e -^(A/)(0)dt 
Jo Jo Jo 

roo 

= / e-^Q 2 s ds-(R p f)(Q). 
Jo 



(6.36) Q^2C 2 +2D 2 sup E(Ma)| 2 ), 

aGR d 

the first part of the lemma follows from the Minkowski- type bound 

(6.37) (|a| + |6|) 1/2 < M^ + I&l 1 / 2 , 

valid for all a, b 6 R. 

The first follows from the second as in the proof of Lemma 5.4. □ 

According to Theorem 1.2, if uq and vo are bounded and measurable 
functions, then there exist a.s.-unique mild solutions u and v to the stochas- 
tic heat equation (6.1). Our next lemma shows that if uq and vq are suitably 
close, then u and v are as well. 

Lemma 6.12. Assume Condition 1.1 holds. Let u and v be two mild 
solutions to (6.1) whose initial functions uo,vo : R d — > R are both bounded 
and in L 1 (R d ). Then there exists /?o which depends only on f and Lip CT; 
such that for all (3 > (3q, 

(6-38) MAu - v) < ^ ^ fj + 2Rem) df. 

Proof. Throughout we suppose that c > 1. Now let us consider a 
fixed x € R d and i > 0, and define u n and to be the respective Picard 
approximation to u and v at the nth stage. Then, 

(6.39) K(z) - <(x)|| 2 < 2 |(P t u )(z) - (P t vo){x)\ 2 + 2QT 1 , 
where Q™ -1 is defined as 

(6.40) f [ p t - s (y - x) (o-(u n -\y)) ~ <v n s -\y))) F{Asdy) 
JO JR d 

Define for all t > and n > 1, 

(6.41) fl-:= supE(|<(a)-<(a)| 2 ). 
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Then it follows from our construction of stochastic convolutions [and (6.19)] 
that Q" _1 is bounded above by 



(6.42) 



u pI / H l l l ds I d y / dz Pt- S (y - x)p t -s{z - x)f(z - 

JO Jn d JR d 



y) 



= Lip 2 - f HrHPt-sfMds. 
Jo 

And hence, 

poo poo 

(6.43) / e-^Q™- 1 dt < Lip 2 • / e-^Hl 1 " 1 dt ■ (Rpf)(0). 
Jo Jo 

This and (6.39) together prove that 

poo poo 

/ e-^H?dt < 2 / sup \(P t u )(x) - (P t vo)(x)\ 2 dt 

(6.44) J ° J ° 

v ' poo 

+ 2Lip 2 • / e-^r'di-^/XO). 
J o 

Condition 1.1 and the monotone convergence theorem together imply that 
lirng_ ) . oo (.R j g/)(0) = 0. We can choose fto so large that 

(6.45) 4Lip2( J R /3o /)(0) < 1. 
Since 

poo 

(6.46) / e-^H? dt = {Mp{u n - v n )} 2 , 
Jo 

it follows that for all ft > fto, 

(poo \ 
/ e-^ sup \(P t u )(x) - (P t v )(x)\ 2 dt 
JO ' H' y 



cGR d 

"W 2 sup \(P t u )(x) - (P t v )(x)\ dt, 
xe~R d 

thanks to Minkowski's inequality. Recall from the proof of Theorem of 1.2 
that, among other things, 



£2 / 

Jo 



(6.48) \\u n - u\\ 2 + \\v n - v || 2 ->■ as n ->■ oo. 
Therefore, Fatou's lemma implies that for all ft > fto, 

poo 

(6.49) N p (u-v)<2 e-^ 2 sup |(P t uo)0c) - (^o)WI dt. 

^0 a;GR d 
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In order to complete the lemma, we choose f3 as large as needed in order 
to ensure that the preceding inequality holds, and define 

(6.50) q :=uq- vq. 

Hawkes's theorem [Proposition 2.3, p. 19] assures us that pt and pt are both 
in L 1 (R d ); in Fourier-analytic language, pt is in the Fourier [or Wiener] 
algebra 

(6.51) A(R d ) := {/i G L 1 (R d ) : h G L 1 (R d )| , 

for all i > 0. It is easy to see that the Fourier transform is one-to-one and 
onto on A(R d ). And continuous elements of A(R d ) are in fact in Co(R rf ), 
by the Riemann-Lebesgue lemma. Since uq and vq are integrable, it follows 
also that PtUo,PtVo G A(R d ) for all t > 0, and that PtUo and PfVo are 
both continuous. The inversion theorem of Fourier analysis can be applied 
pointwise to the elements of ^4(R d ). Therefore, 

sup \(P t q)(x)\ < ||^|| L i (Rd) 

(6.52) =<£?U*-"™*®\* 



*(2 



tt) Jr* 



And hence, 



(6-53) ire-^su P \(P tq )( X )\ d t<'[ '^f } ~ de 
2 y ieRl! (2vr) (1 y Rd /3 + 2ReW(f ) 

This proves the lemma. □ 



6.3. Proof of Existence 

Having prepared the background material, we are ready to demonstrate 
Theorem 6.2. Before we begin the proof, let us recall the following: 

Definition 6.13. Let {ip n }^ =1 denote a collection of measurable func- 
tions from R rf to R + such that: (i) n h-> ||Y'n||Li(Rd) is uniformly bounded; 
and (ii) lim n _ yoo ip n = I pointwise. Then we say that {ip n }^ =1 is a weak 
mollifier. 

Clearly, weak mollifiers are mollifiers in the usual sense. But the converse 
is not in general true, because weak mollifiers need not have very good 
smoothness properties. 
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We need this definition because our strategy of the proof is the following: 
Let {V'nj^Li denote a weak mollifier. Then, according to Theorem 1.2, we 
can solve (6.1) subject to the initial condition uq * ip n , since uq * ip n is a 
bounded and measurable function. If denotes the solution, we then 
proceed to show that u := linin^oo exists in a suitable sense and is a 
temperate solution to (6.1). Now we write down the details of this argument. 



PROOF of Theorem 6.2. Recall that uq is a finite Borel measure and 
define := ip n * uo, where {^nj^Li is a weak mollifier. According to 
Theorem 1.2, there is a solution to the following SPDE: 

(6.54) = 0&4 B) )(*) + ^ui n \x))F t (x), 

and the solution is unique up to evanescence. Theorem 1.2 also assures us 
that 

(6.55) C a := sup sup e~ at E ( (x) ) < oo, 
provided that a > is sufficiently large. Note that if 6 > a, then 

poo 

Ne{vS n) )<C a e-V-^dt 

(6.56) h 

= < oo. 

V — a 

That is, £ Lg for all 9 sufficiently large. According to Lemma 6.12, for 
all n,m> 1, 

(6.57) N, (.(.. - .0-)) < ^ J^L . \ U Q _ d { , 

provided that /3 > (3q for a /3q that depends only on / and Lip CT . Therefore, 
it follows that if (3 > (3q, then {u^}^ =1 is a Cauchy sequence in L^. Let 
u := {ut(%)}t>o,xeK d denote the limit. By definition, 

(6.58) lim «(") = u in Pi Ll. 

/3>/3o 

And therefore, Lemma 6.11 and our extension of stochastic convolutions 
together imply that 

(6.59) lim ( - Au^ n A =u-Au in Pi Ll, 
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where Au^ is defined in (5.5) on page 56; and Au is denned in the same 
way, but now interpreted as a stochastic convolution in the sense of the 
present chapter. 

In particular, for all T > 0, 

(6.60) [ T sup E (\[u r t \x) - (Au n ) t (x)} - [u t (x) - (Au) t (x)}\ 2 ) dt -»• 0, 

Jo x£R d V J 

as n — > oo. This follows simply because 

rT /*oo 

(6.61) / K(s)ds<e^ T e- /3s K(s)ds, 
Jo Jo 

for all nonnegative measurable k : R+ — > R+. 

Since is a mild solution to (6.1) with initial data = ip n * no, 

Tonelli's theorem implies that 

(6.62) ut\x) - (Au n ) t (x) = ((P t u ) * ^n) (x). 

Hawkes's theorem [Proposition 2.3, p. 19] implies that P t UQ is uniformly 
continuous for every t > 0. Therefore, for every t > fixed, 

(6.63) lim (-PfUo) * V'n = -ft^o uniformly. 

n— >oo 

It follows from the preceding and Fatou's lemma that for all T > 0, 

(6.64) f sup E (\u t (x) - (P t u )(x) - (Au) t (x)\ 2 ) dt = 0. 
Jo xe~R d v ' 

This proves the theorem. □ 

The following is the final result of this section. It shows that the temper- 
ate solution to (6.1) is unique among a natural family of possible solutions. 
It is entirely possible that one can establish the uniqueness of the temperate 
solution among a larger family than that offered below. But we are not 
aware of such further improvements. 

PROPOSITION 6.14. The temperate solution u, provided by the preceding 
proof, does not depend on the choice of the weak mollifier, up to evanescence. 

Proof. Let {K n }'^ ) = i be another weak mollifier, and define v to be the 
temperate solution that the preceding proof yields based on {k„}^ =1 ; that 
is, v n is obtained from {k/} / c ^ 1 in the same way that u n was obtained from 

Let (3o be as in the proof of Theorem 6.2 [see also Lemma 6.12], and 
recall that fto does not depend on {ip n }^=i, {^n}^Li, etc. In accord with 
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Fatou's lemma, 

Mr(u — v) < liminf Mr {u 

n— >oo 

(6 - 65) < 4 /• |n (OI 



liminf / '^ (0 ' • - M0 d£, 



(27r) d 

for every /3 > po, owing to Lemma 6.12. And the preceding is zero by the 
dominated convergence theorem. This proves the proposition. □ 



6.4. More A Priori Estimates 

Let uq be a finite Borel measure on R rf that satisfies (6.10). In order 
to investigate the large-time behavior of the temperate solution to (6.1), we 
need to introduce a suitable family of Banach spaces. The task of the present 
section is precisely to do that. Our analysis of the temperate solution to 
(6.1) will resume after this section. 

Define for all j3 > 0, finite real numbers p > 1, and predictable random 
fields v := {v t (x)} t>0jXen d, 



(6.66) 




so that Npp (v) is the same quantity as Mr (v) ; the latter was defined earlier 
in (6.21). And, just as was the case with Mr when p = 2, every Mp tP is a 
norm on equivalence classes of p-times integrable predictable processes that 
are modifications of one another. 

Definition 6.15. Let Lip denote the collection of all predictable random 
fields V := {^t(^)}f>0,a;GR d sucn tnat Nf3 )P (v) < oo. 

Thus, our present definition of L 2 ^ agrees with our older one. 
The following generalizes Lemma 6.11 to the case that p is an integer 
> 2. 

Lemma 6.16. For all j3 > 0, even integers p > 2, and Walsh-integrable 
random fields v and w, 

(6.67) Mr, p (Av) < z v + D CT A/»)) 7^7X0), 
and 

(6.68) Mp, p (Av - Aw) < z p U Va Mp, p {v - w)yj2{R p f)(Q). 
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Proof. The ensuing argument is a generalization of the proof of Lemma 
5.4. Indeed, by (5.20), E(\(Av) t (x)\ p ) is at most z v p times 

[ ds [ dy j dz V^(y , z)p t ^ s (y - x)p t - s (z - x)f(y - z)) 
Jo Jn d Jn d J 



(6.69) 
where 



V;(y,z) := \\a (v s (y))\\ p ■ \\a (v s (z))\\ p 



(6-70) < [C a + B a sup \\v s (q)\\ p 

V q€K d . 



<2G S , 



where 



(6.71) Q s :=C 2 a + Bl sup \\v s (q)\\ 2 p . 

It follows that E(\(Av) t (x)\ p ) is bounded above by 2 p / 2 zf, times 

E(\(Av) t ( X )n 

(6 72) < 2 p/2 ^ ^ Qs ds dy dz p t _ s {y)p t _ s (z)f{y - z)j 

/ft \ P/2 

= 2 p/2 ^fy o e s (p t _ s /)(o)d s J , 

thanks a direct computation that was made already during the course of the 
proof of Lemma 5.4. 

Let us examine the preceding display next: The right-most quantity is 
independent of the variable x. Therefore, we can maximize the extreme 
terms in that display over all x G R rf , raise the resulting inequality to the 
power 2/p on all sides, and then finally multiply by exp(— fit) and integrate 
[dt] to obtain the following: 



(6.73) Mj3,p{v) < %^2^ e-^0 s ds • (^/)(0) 

The first bound of the lemma follows after a few direct computations. The 
second bound follows from the first in a manner that has been pointed out 
several times already; see, for example, the end of the proof of Lemma 5.4 
for an outline. □ 

Note that if v is a Walsh-integrable random field, then so is aov, and then 
Av is none other than the stochastic convolution p * (a o v)F. Therefore, 



6.4. MORE A PRIORI ESTIMATES 



91 



the following is a refinement of Lemma 6.16 that works in the setting of 
stochastic convolutions. 

Lemma 6.17. Recall that a : R — > R is a predetermined nonrandom and 
Lipschitz- continuous function. Choose and fix a real (3 > and an even 
integer p > 2. Then, for all predictable random fields v and w, 

(6.74) N/3,p [p * (a o v)F) < z p + D CT A/»)) y/WpfM, 
and 

Mp p (p * (cr o v)F — p * (a o w)F J 

(6.75) ' V J 

< z p Up a Mf3, p (v - w)sj2{R p f){Q). 

PROOF. We can, and will, assume without loss of generality that the 
three quantities (Rpf)(0), Afp :P (v), and Afp jP (v — w) are all finite. 

Since p > 2, it follows that the stochastic convolution p * vF is well 
defined. Moreover, we can find a sequence w 1 ,^ 2 ,... of Walsh-integrable 
random fields such that p * v n F converges to p * vF in L». In particular, 
there exists a Lebesgue-null set Nq C R+ such that for all t No, 

(6.76) lim sup E (\v?{x) - v t {x)\ 2 ) = 0. 

Therefore, we may apply Fatou's lemma to the [pseudo-] norm ftfp ;P and 
deduce the following: 

Mg v \p* vF) < lim inf Ms v (p* v n F) 

(6.77) ' V 7 n ^°° ' V ' 

< ^lirninfA^, p (t; n ) y/(R p f)(0), 

because we can apply Lemma 6.16 to each Walsh-integrable process v n [with 
a(u) := u]. Among other things, this proves that p * vF G L^, and 

(6.78) A/>, p (p * vF) < Zp Mp, p (v)^2(Rpf)(0). 

The first inequality of the lemma follows from the above and Minkowski's 
inequality, since a(vt(x)) < C a + D a \vt(x)\ pointwise. 

Because we also have that Np :P (w) < oo, we can deduce the second 
inequality from Lemma 6.16 as well. □ 



We will not need the entire strength of the preceding lemma itself. We 
have stated it for two reasons: First of all, it suggests that our extension 
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of the stochastic convolution has good continuity properties, viewed as ele- 
ments of the Banach spaces L p ^; and also, this is a simple setting in which 
an approximation method introduced via the proof of Lemma 6.17. We will 
use that method later on. 



6.5. An Upper Bound for Growth 

The goal of this section is to establish Theorem 6.7. 

Proof of Theorem 6.7. Let {ifj n }^ =1 be a weak mollifier, and let 
u {n) ._ {u[ n \x)} t>0 xg Rd be the solution to (6.1) with initial data ip n * uq. 
By the Minkowski inequality, the following is valid for every integer n > 1: 

(6.79) A/>, p (ttW) < Mp, v (P.uo) + Mfs, p (p * (a o U W) F) . 

The first term on the right is estimated easily, thanks to Minkowski's in- 
equality, as follows: 

/ oo \ l l 2 

A/>, p (P.u ) = ( fe"^ sup \(P t u )(x)\ 2 dt) 
\Jo xeR d J 

( 6 - 8 °) < He-^ 2 sup \{P t u )(x)\ dt 

Jo xen d 



" (2vr) d L 



x£~R d 

\M0\ 



l Rd p + 2Re\P(£) 
see (6.53). Therefore, Lemma 6.17 implies that 



(6.81) 



\f (,An)\ < 2 f \M0\ j t 



+ z P (jfc + D^, p («(»>)) y/2(Rpf)(0). 



Let us assume, for the time being, that we could prove that Afp p (u^) is 
finite — that is, is in the Banach space L p ^ — provided that 2 

(6.82) z^D a ^2{R p f){Q) < 1. 

Then we could rearrange (6.81) and find that 

(6 , 3) ^^^l-J^v + ffl, 



'(6.82) holds for all j3 sufficiently large because thanks to (1.1), fim / g_ >oo (.R / g/)(0) = 0. 



6.5. AN UPPER BOUND FOR GROWTH 



9.3 



where 

(6.84) c := 1 

l-ZpDefaRpfW) 

And therefore, the proof of Lemma 6.17 implies that 
Mfa{v) <liminfA^, p ( W W) 

n— i-oo 

(6 - 85> <J*f !Mli df + aM 



In other words, we have shown the following: If we could prove that € 
L P p for every n, then in fact {■u^ n ^}^L 1 is bounded in Ll whence u £ L p ^ 
and Theorem 6.7 follows. 

Now we proceed as follows: Let 

(6.86) ^ n,0) (x) := (ip n * u )(x) for all t > and z € R d , 
and iteratively define 

(6.87) u^ k \x) := (P t u )(x) + J* J p t _ s (y-x)a (u^'*" 1 ^)) F(dyds). 
By the Minkowski inequality and (6.53), 

A/>, p (u^°A = ( fV* sup |(P t (V re *«o)) (x)| 2 dt) 

(6.88) / 

< 2 /• |^(0|-|^o(0| 

" (2vr^y Rd /3 + 2Re*(0 ? ' 

which is finite; it is in fact bounded uniformly in n since sup n>1 ||V'n||L 1 (R d ) 
is finite by the very definition of weak mollifiers. Now the argument that 
led to (6.81) leads also to the following: 

J\f R ( u (n,k+lA < _J_ f M^j d £ 

(6g9) N "' P[ ~ J-(2^y Rd /? + 2Re^(O de 



+ z % 



(w~2 + D ^p { u(n ' k) )) y/^Rpfm- 



V/3V2 

This, (6.88), and induction together prove that if (6.82) is in effect, then 

(6.90) supA/>, p (V n ' fe) ) < oo. 

k>i ' ^ ' 

By the Borel-Cantelli lemma, and owing to Theorem 1.2 and its proof, 

(6.91) lim sup E ( up k \x) - uj n) (x) ) = for all t > 0. 
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[In fact, this is essentially (6.76), but we have noted further that "almost 
all t" can be replaced by "all i" in the present setting, since is a mild 
solution to (6.1).] Therefore, Fatou's lemma proves that 

(6.92) Mb p (u {n) ) < liminf Mr p (u^ k A < oo. 

This establishes Theorem 6.7. □ 
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